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I.  INTRODUCTION 

Considerable  interest  has  been  shown  in  the  study  of  residence-times  in 
chemical  equipment.  This  is  due  to  the  fact  that  the  residence-time  distri- 
bution gives  valuable  information  on  the  nature  of  the  flow  through  a  plant. 
The  total  quantity  of  material  inside  a  working  plant  is  often  difficult  to 
determine  accurately,  and  a  calculation  of  the  mean  residence-time  is  valuable 
in  giving  the  material  inventory  of  a  plant. 

Consider  the  case  of  chemical  reactions  which,  if  not  zero  order  or  auto- 
catalytic,  decrease  in  rate  as  time  proceeds.  It  is  possible  to  calculate 
the  average  yield  of  the  reaction  by  means  of  the  following  equation  (1) 

-2_=J^f(8)E(6)dO  (1-1) 

C0    0 

where  C  is  the  mean  concentration  of  reagent  left  in  the  effluent,  Cn  is  the 
concentration  of  reagent  in  the  feed,  f(9)  is  the  kinetic  equation  which  gives 
C/CQ  as  a  function  of  time  in  a  batch  reaction,  and  E(6)  is  the  residence-time 
distribution  function.  This  expression  is  only  correct  for  first  order 
chemical  reactions.  However,  the  residence-time  study  provides  valuable 
information  for  the  explanation  of  reactor  behavior  even  in  cases  where  the 
reaction  is  not  first  order. 

Another  field  of  fluid  flow  which  is  increasing  in  importance  industrially 
is  the  behavior  of  fluids  which  do  not  obey  Newton's  law  of  viscosity  in 
motion.  Because  of  this  negative  definition  of  non-Newtonian  behavior,  no 
single  equation  can  describe  exactly  the  shear-stress  and  shear-rate  relation- 
ships of  all  such  materials  over  all  ranges  of  shear  rates.  Numerous  empirical 
equations  have  been  proposed  to  express  the  steady  state  relation  between  the 
shear- stress  and  shear- rate.  Among  these  the  Ostwald-de  Waele  model  and  the 


Bingham  plastic  model  are  the  most  generally  useful  two-constant  models,  and 
the  Ellis  model  is  the  simplest  three-constant  model. 

Most  of  the  investigations  of  fluid  dispersion  have  been  confined  to 
Newtonian  fluids.  But  fluid  dispersion  in  the  processing  of  non-Newtonian 
fluids  is  as  important  as  it  is  in  the  processing  of  Newtonian  fluids.  This 
thesis  will  be  concerned  mainly  with  the  development  of  mathematical  models 
which  characterize  the  dispersion  of  non-Newtonian  fluids  in  flow  systems. 

A  generalized  mathematical  expression  for  diffusion  in  a  flow  system  at 
constant  temperature  and  pressure  is  (2) 

!£-+  V-VG  =DV2Cf  R  (1-2) 

where  D  is  assumed  to  be  independent  of  C. 

Suppose  that  a  fluid  flow  through  a  cylindrical  tube  with  the  symmetrical 
concentration  distribution  about  the  central  line  of  the  tube.  When  there  is 
no  chemical  reaction  in  the  system,  Equation  (2)  becomes 

S-^+iS+S^xWg  (1"3) 

It  would  be  difficult  to  find  a  complete  solution  of  Equation  (3)  which 
gives  the  value  of  C  for  all  values  of  r,  x  and  t  when  the  distribution  C  at 
time  t=0  is  known.  However,  approximate  solutions  can  be  found  which  are 
valid  in  certain  limiting  conditions.  When  the  effects  of  both  cross-sectional 
and  longitudinal  diffusions  are  negligible,  the  steady-state  velocity  profile 
becomes  the  only  factor  governing  the  apparent  fluid  dispersion.  The  flow 
model  corresponding  to  such  a  condition  is  usually  called  the  convective  model. 
When  the  variation  of  axial  velocity  with  cross-sectional  position  and  the 
cross- sectional  material  transport  by  molecular  diffusion  are  assumed  to  be 
the  dominant  dispersion  mechanisms,  the  solution  of  Equation  (3)  is  called 


the  dispersion  model.  Both  models  are  treated  in  this  thesis  for  various 
non-Newtonian  velocity  distribution  functions,  such  as  Bingham  plastic, 
Ostwald-de  Waele  and  Ellis  models  and  various  flow  geometries,  such  as 
cylindrical  tubes,  slits  and  open  channels. 

The  analysis  of  turbulent  flow  of  non-Newtonian  fluids  through  smooth 
round  tubes  was  performed  for  the  first  time  by  Dodge  and  Metzner  (3).  The 
analysis  permitted  the  prediction  of  non-Newtonian  turbulent  velocity  pro- 
files. The  dispersion  model  of  turbulent  non-Newtonian  flow  is  also  dis- 
cussed in  this  thesis. 


II.  LITERATURE  SURVEY 

Fluid  is  passed  through  process  equipment  so  that  it  may  be  modified  in 
one  way  or  another.  To  predict  the  performance  of  equipment,  we  must  know 
the  complete  flow  pattern  of  the  fluid  within  the  vessel.  Because  of  the 
practical  difficulties  connected  with  obtaining  and  interpreting  such  infor- 
mation, an  alternate  approach  is  used,  which  requires  knowledge  only  of  how 
long  different  elements  of  fluid  remain  in  the  vessel.  Though  this  partial 
information  is  not  sufficient  to  completely  define  the  nonideal  flow  within 
the  vessel,  it  is  relatively  simple  to  obtain.  The  field  today  can  be 
divided  into  several  areas.  One  of  these  is  how  residence-time  distribution 
functions  can  be  measured  experimentally.  Another  area  which  has  had  consid- 
erable attention  is  the  theoretical  derivation  of  the  relationships  between 
effective  axial  dispersion  coefficients  and  dispersion  models.  It  is  the 
purpose  of  this  chapter  to  review  the  present  situation  in  this  field. 

The  convective  model  of  dispersion  for  Nex-rtonian  fluids  in  laminar  flow 
was  proposed  by  many  investigators  (4-,  5)  to  predict  residence-time  distribu- 
tion functions  for  flow  in  circular  conduits.  The  cummulative  age  distribu- 
tion at  the  outlet  of  the  system, ;  which  corresponds  to  the  response  to  a 
step  function  input  of  a  tracer  without  entrance  effect  being  considered,  is 
expressed  as 


F(e)=  i  — L-  for  e^i 

(20)2  2 


=  0  for  e<i 

""  2 


(2-1) 


where  Q   is  the  dimensionless  time  defined  by  0  =  tVx/L.  The  exit  age  distri- 
bution, which  corresponds  to  the  response  to  a  Dirac  delta  function  input  of 


a  tracer,  is 


1(9)=  1  for  e^- 

■     I?  2 


=  0  for  0^i 

2 


(2-2) 


If,  however,  the  entrance  effect  cannot  be  neglected,  the  initial 
distribution  of  the  tracer  is  uniform  spacewise  along  the  tube.  For  this 
case  the  curaraulative  and  exit  age  distributions  were  found  to  be 


F(9)=  1  -  JL  for  e^i 

-W       2  8  2 


=  0  for  0<  I 

2 

S(9)=  -i-  for  Q>± 

ze^  2 


=  0  for  0^- 

2 


(2-3) 


(2-^) 


In  order  to  give  a  more  precise  description  of  the  dispersion  character- 
istics of  fluids,  axial  and  radial  diffusion  and  the  velocity  profile  should 
be  considered  simultaneously.  In  a  series  of  papers  (6,  7,  8),  Taylor  has 
treated  the  dispersion  of  soluble  natter  in  a  solvent  flowing  through  a 
circular  tube.  For  a  Newtonian  fluid  in  laminar  flow,  the  distribution  of 
concentration,  C,  of  the  soluble  material  depends  on  the  balance  between 
convection  along  the  tube  due  to  variation  in  velocity  over  the  cross  section 
and  on  cross-sectional  molecular  diffusion.  The  corresponding  partial 
differential  equation  is 

C(S  +  S  +  r  5?)=  S?  +  V»(1  "  i]  M  (2-5) 


Here,  D,  the  coefficient  of  molecular  diffusion,  is  assumed  to  be  independent 
of  C.  In  general,  the  transfer  of  C  along  the  tube  by  molecular  diffusion  is 

small  compared  with  that  produced  by  convection.  It  is  thus  assumed  that 

2  2 

—  C  is  negligible  compared  with  -2_ Q.  a.  1  _?_£   The  condition  necessary  for 
ax2  ar2   r  Sr. 

this  to  be  true  can  be  expressed  as  (8) 

R2?2 
D  >>  £_  (2_6) 

48D 
The  transport  equation  therefore  takes  the  form 

subject  to  the  boundary  conditions  that  the  wall  of  the  tube  is  impermeable 
and  that  the  concentration  distribution  is  symmetrical  with  respect  to  the 
center  line  of  the  tube,  or 

JL£=0,  at  r  =  R  and  r  =  0    (2-8) 

5  r 

It  would  be  difficult  to  find  the  complete  solution  of  Equation  (7); 
on  the  other  hand,  certain  approximate  solutions  do  exist  which  are  valid 
under  certain  limiting  conditions.  One  such  approximate  solution  considers 
dispersion  by  convection  alone,  which  applies  to  the  case  when  both  the  cross- 
sectional  diffusion  and  longitudinal  diffusion  are  negligible,  and  the  steady 
state  velocity  profile  becomes  the  only  factor  governing  the  apparent  fluid 
dispersion.  Some  calculated  distributions  of  G  along  a  tube  have  been  given 
by  Taylor,  which  can  be  reduced  to  the  residence-time  distribution  functions 
given  in  Equations  (1)  through  (4).  Another  approximate  solution  suggested 
by  Taylor  is  that  the  time  necessary  for  a  radial  variation  in  C  to  die  down 
owing  to  radial  diffusion  is  much  shorter  than  the  time  necessary  for  an. 
appreciable  change  in  C  to  occur  through  longitudinal  convection.  This  can 


be  expressed  by  the  condition  (6) 

_  >>    2r2  (2-9) 

Vx      (3.8)2D 

where  V  is  the  mean  speed  of  flow  and  L  is  the  longitudinal  extent  of  the 
x 

region  in  which  dC/3x  is  appreciable. 

It  is  convenient  in  the  present  discussion  to  define  concentration  and 
velocity  relative  to. axes  which  move  with  the  mean  flow.  The  velocity 
relative  to  these  axes  is 

V  =  2V  (1  -  £-)  -  V  =  V  (1  -  2E- )  (2-10) 

xi     x'    R2    x    xs     R2 

Letting  ^  =  r/R,  the  equation  of  diffusion  becomes 

2- 
32C  ,    1    SC_R2   ^C    ,     '    x/n        o*2nSC  (9  T-.N 

"^2+5^"F7t  +  ~u  "  *  ;Tx~  K    1±) 

where  x,  =  x  -Vxt.  For  small  radial  variations  in  C,  partial  equilibrium  may 

be  assumed,  that  is,  the  rate  of  change  with  time  is  equated  to  zero.  Since 

the  mean  velocity  across  planes  for  which  X-,  is  constant  is  zero,  the  transfer 

of  C  across  such  planes  depends  only  on  the  radial  variation  of  C.  In  this 

calculation,  ^C/^x.  is  taken  to  be  independent  of  £.  The  small  radial 

variation  in  C  can  therefore  be  calculated  from  the  equation 

2- 

±!c  +  iac  =  i\(1  „  2|2)  |_c_  (2_12) 

^|2  1*1         V  ^xx 

Since  ^C^x.,  is  assumed  to  be  independent  of  £,  and  if  C  is  the  mean 
concentration  over  a  section,  3 c/dx,  is  indistinguishable  from  ^Cm/&Xn . 
Therefore,  Equation  (12)  can  be  solved  to  give 
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The  rate  at  which  C  is  transported  across  a  section  at  x,  is 


Q=  2TfR2J*  Vx(l  -  2^2)C^d£  (2-14) 


-1 
0 
Inserting  the  value  of  C  from  Equation  (13).  Equation  (14)  is  found  to  be 


2  R2^2  3fi, 

Taylor  has  concluded  that  the  combined  effects  of  longitudinal  convection 
and  radial  molecular  diffusion  are  equivalent  to  the  transfer  of  the  solute 
across  planes  which  move  with  the  mean  speed  of  the  flow  and  the  rate  of  this 
transfer  is  equal  to  that  which  diffusivity,  K,  would  give  to  a  stationary 
fluid,  as  shown  by 


Tr2 


K=  £  (2-16) 

48D 

The  analogous  problem  of  dispersion  in  turbulent  flow  can  be  solved  in 
the  same  manner  (7).  Taylor  has  assumed  in  his  analysis  that  the  universal 
velocity  distribution  in  a  pipe  and  Reynolds '  analogy  that  the  transfer  of 
matter,  heat,  and  momentum  by  turbulence  are  exactly  analogous.  The  apparent 
diffusion  coefficient  K  under  this  case  has  been  found  to  be  10.1RV\  or 
K=  7.14RV  If.     Here  R  is  the  radius  of  the  pipe,  Vx   is  the  mean  flow  velocity, 
f  is  the  Fanning  friction  factor  and  V^  the  friction  velocity. 

In  a  later  paper  {9)   Aris  has  presented  a  new  analysis  in  which  the 
restrictions  imposed  by  Taylor  are  removed.  He  obtained  this  by  fixing 
attention  on  the  movement  of  the  center  of  gravity  of  the  distribution  of 
solute  and  the  growth  of  its  higher  moments.  It  has  been  shown  that  the  rate 

of  growth  of  the  variance  is  proportional  to  the  sum  of  the  molecular  diffusion 

2—2 

coefficient  D,  and  the  Taylor  dispersion  coefficient  KR  V  /D.  The  value  of  K 

for  the  Newtonian  system  has  been  given  by  Taylor  to  be  l/48  for  laminar  flow 


and  10.1  for  turbulent  flow  (6,  7). 

Aris  (10)  has  also  applied  this  method  to  the  case  in  which  the  solute 
can  pass  into  another  fluid  phase  flowing  through  a  tube.  The  apparent 
diffusion  coefficient  has  been  shown  to  be  the  sun  of  the  molecular  diffusion 
coefficient  and  the  Taylor  dispersion  coefficients  in  the  t;/o  phases,  and  a 
term  due  to  the  finite  rate  of  partition  between  them.  In  another  paper  (11), 
Aris  has  treated  the  case  of  a  viscous  flow  under  a  pulsating  pressure 
gradient.  It  has  been  found  that  the  Taylor  diffusion  coefficient  contains 
terns  proportional  to  the  square  of  the  amplitude  of  the  pressure  pulsations. 
But  the  coefficients  of  such  terms  rarely  contribute  more  than  a  fraction  of 
l/l28  to  the  total  dispersion  coefficient. 

Taylor's  treatment  of  turbulent  flow  is  valid  only  for  high  Reynolds 
numbers,  however,  because  the  velocity  profile  used  in  the  treatment  is  valid 
only  when  the  laminar  sublayer  and  transition  layers  are  negligibly  small. 
Tichacek  and  co-workers  (12)  have  refined  Taylor's  method  by  including  the 
effect  of  molecular  diffusion  and  by  introducing  the  experimental  velocity 
profiles  rather  than  a  generalized  profile.  They  have  included  a  first  order 
approximation  to  9C/dxas  a  function  of  the  profile.  They  have  also  shown 
that  the  effects  of  axial  turbulent  diffusion  are  negligible  compared  to  the 
mixing  caused  by  radial  differences  in  the  velocity.  According  to  their 
analysis,  Dj/V  d  is  dependent  of  the  Reynolds  number,  pipe  roughness,  and 
Schmidt  number.  They  have  reported  that  their  theoretical  data  are  applicable 
with  less  than  25$  error  as  shown  by  comparison  with  Taylor's  experimental 
data. 

Giddings  and  Seager  (13)  have  theoretically  and  experimentally  investi- 
gated a  method  siiailar  in  operation  to  chromatography  techniques  for  measuring 
a  wide  range  of  diffusion  coefficients.  Their  experimental  work  deals  with 
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gaseous  diffusion  coefficients  measured  at  various  flow  velocities,  concen- 
trations, etc.  The  final  equation  is 

where  h  is  the  height  equivalent  to  a  theoretical  plate  in  a  typical  gas 
chromatograph  and  the  right-hand  side  is  the  sum  of  the  molecular  diffusion 
coefficient  and  Taylor's  dispersion  coefficient. 

Bailey  and  Gogarty  (1^)  have  performed  experiments  in  which  a  dilute 
solution  of  potassium  permanganate  is  displaced  by  water.  They  also  developed 
a  numerical  method  and  obtained  an  accurate  solution  of  the  mixing  equation 
when  longitudinal  diffusion  is  neglected.  They  have  reported  that  good  agree- 
ment does  exist  between  the  numerical  solution  and  the  given  experimental 
results,  that  the  mixing  zone  lengths  resulting  from  their  experiment  are 
nearly  the  same  as  those  predicted  by  Taylor  for  a  narrow  range  of  dimension- 
less  time,  6  =  Dt/R  .  However,  experiments  at  various  flow  times  with  a 
fixed  velocity  have  shown  that  experimental  dispersion  coefficients  increased 
slowly  with  time. 

Bournia  and  co-workers  (15)  have  measured  the  longitudinal  dispersion  of 
a  finite  slug  of  gas  at  various  velocities  by  using  a  gas  (1,3-butadiene)  that 
absorbs  light  in  the  ultra-violet  region  and  passing  the  dispersed  slug  through 
a  narrow  beam  of  ultra-violet  light  of  wave-length  250  m/A.  The  results 
indicate  that  Taylor's  approximation  does  not  apply  for  values  of  Vx  below 
2  cm/sec,  and  that  the  criterion  O/D  »  6.9  is  very  important  while  the 
upper  criterion  4L/R  »  VXR/D  is  less  important  because  Taylor's  approximation 
appears  to  apply  best  at  the  higher  velocities. 

Farrell  and  Leonard  (16)  also  derived  a  solution  for  the  laminar  flow 
problem  when  the  axial  molecular  diffusion  is  neglected.  Their  solution 
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involved  a  series  of  cigen  functions  of  the  Laplace  transformation  of  the 

concentration.  The  moments  of  the  concentration  distribution  can  readily  be 
evaluated  by  making  use  of  a  computer. 

Taylor  and  Aris  both  have  concluded  that  an  effective  axial  dispersion 
coefficient  K  can  also  be  used  in  turbulent  flow.  This  coefficient  has  been 
found  to  be  a  function  of  the  well-known  Fanning  friction  factor.  Bischoff 
and  Levenspiel  (17)  have  extended  Aris'  theory  to  include  a  linear  rate 
process,  and  have  used  the  result  to  construct  comprehensive  correlations  of 
dispersion  coefficients. 

Hawthorn  (18)  has  considered  the  temperature  effect  of  viscosity  on  the 
dispersion  coefficient  and  has  found  that  it  can  be  altered  by  a  factor  of 
two  in  laminar  flow,  but  that  there  is  little  effect  for  fully  developed 
turbulent  flow. 

The  analysis  of  dispersion  accompanying  the  flow  of  a  non-Newtonian  fluid 
has  recently  been  considered  by  Fan  and  Hwang  (19).  They  have  derived  convec- 
tive  models  for  Bingham  plastic  and  Ostwald-de  Waele  fluids.  The  dispersion 
model  for  Ostwald-de  Waele  fluids  flowing  through  a  circular  tube  has  also 
been  derived  for  which  Taylor's  results  can  be  considered  as  a  special  case  of 
Fan  and  Hwang's  expressions. 
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III.  CONVECTIVE  MODELS 

A  general  discussion  of  residence-tine  distributions  for  some  flow 
models  has  been  made  by  Fan  and  Hwang  (19).  First  of  all,  they  have  treated 
the  convective  model  extensively  for  fluids  which  can  be  represented  by  the 
Ostwald-de  Waele  and  Bingham  plastic  models.  In  this  chapter,  the  similarity 
of  residence-time  distributions  between  fluid  flow  through  slits  (opening 
between  two  parallel  plates)  and  fluid  flow  through  open  channels  iri.ll  be 
verified.  The  convective  model  for  the  laminar  flow  of  the  Bingham  plastic 
fluid  which  is  more  general  than  that  derived  by  Fan  and  Hwang  (19)  will  also 
be  derived.  Constant  molecular  diffusion  coefficients,  steady  and  iso- thermal 
flow  situations,  and  constant  vessel  geometries  will  be  assumed  in  the  present 
analysis. 

FLOW  THROUGH  OPEN  CHANNELS 
1.  Bingham  Plastic  Model 

The  steady- state  rheological  behavior  of  Bingham  plastic  fluids  can  be 
expressed  as  (2) 

€=-{w ~° }t  fori  ($:t)±  T*  (3-D 

A  =  0  fori("C:r)<r  (3-2) 

2  0 

where  >*  and  X  are  parameters  which  characterize  a  fluid. 

Under  conditions  for  which  a  one-dimensional  rheological  statement  in 
rectangular  coordinates  is  valid  (see  Fig.  (l)),  Equations  (1)  and  (2)  reduce 
to 

tyx=T0->\)  J  .  tor    Xyx    >T0  (3-3) 
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dV 

-r-2  =  0  for  -c   <  TA  (3-^) 

dy  yx  v  0 

respectively,  and  the  equation  of  motion  reduces  to  (2) 

JLr  =.pg   cosf  (3-5) 

Mailing  use  of  the  boundary  conditions  at  the  liquid-gas  interface 

X^  -  0  at  y  =  0  (3-6) 

The  equation  of  motion  may  be  integrated  to  give 

"C^-fgy  oosf      .  (3-7) 

Combining  Equations  (3)  and  (7),  and  integrating  the  resulting  equation 
subject  to  the  following  boundary  condition 

Vx=  0  at  y  -S  (3-8) 

one  obtains 

vx  =  -'cs4FLi  +  lr  y  +  ci  C3-9) 

and 

c  =  eEii^os^  _Tis  (3_10) 

1     2^o     ^0 

Hence,  the  velocity  distribution  is 

Vx-  PES*  cos*  [l  -  (g)2]  -Zgf  (1  -  |)  fory>yQ  (3-11) 

0  'O 

m  2^0  l  S       J       h0  S  ~     0 

where 


0       fg  cos^f 
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Introducing  the  dimensionless  distance  variable  £,    =  —  -,   the  velocity  distri- 


bution  nay  also  be  written  as 


Z  ..  .       o 


Vn  =  pcr^  cos^    (1  -  40r  for    y  <  y0  (3-13) 

m  ou  u  0 

^r0 


V    = 


££^2-1   ((l-40)2-(|.40)2] 


2HQ 


,y      >    ^ 


=  Vn  [  1  -  Is -li-  ]  for    y  >  y  (3-14- ) 

d-4o> 

The  mean  velocity  of  the  fluid  is  obtained  by  summing  all  the  velocities 
over  a  cross-section  and  then  dividing  by  the  cross-sectional  area,  that  is, 

J    J    Vx(y)dydz 
K-ffi 0-15) 

J     J       dydz 
0     0 

where  W  is  the  width  of  the  channel.     Therefore 

5l0       m         Jyo»l         (l-i0)2J 
.fK^owy  (j-t0t|{03)  (3-16) 

Dividing  Equation  (16)  by  Equation  (13)  yields 

Vx=|vmp  (3-17) 


wnere 


p  =  i-iio  +  l*o3 . ilia  (3_18) 


•v2 
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When  the  entrance  effect  is  negligible,  the  cumulative  age  distribution, 
■which  is  a  response  to  a  unit  step  function  input,  is 

^e>=^r/'  f   vjy)*&*  (3-19) 


■yo„    ,_       i     ry    „    f  ,      <?-*o> 


2 


=  4-  /  °:ym  dy  +  4-  7      vJi-i^iiLjdy  O-20) 

"x  0         SV*%     '     <l-i„> 
=  ^(9)  +  5^(6)  (3-21) 

where 

F  (9)  =  -1-  J7°  V m  dy  =  ^0  .  Us  (0  -  |p)  (3-22) 

From  Equation  (1*0  one  has 


-V^-Mi^ 


and 

<L  (I)  =   ^(1"^L.  (3-2*+) 


de  * 


Ji- 


2  A  2£ 

36 


Substituting  these  equations  into  Equation  (21)  yields 

_r6   pft-tp) 


#9)«j[ 


ae 


|p  3S3  JTTi 

~    36 


1? 


=  (i  -  u 


for  9  >  r' 


(3-25) 


The  response  to  a  unit  step  function  input  is  therefore 


34, 


f(6)  =  ^.  os(e-s»  +  (i-40) 


, ,  ,f  JH8.  jh1 


3/2       +  I/? 


for    e  >    |p 


(3-26) 


=   0 


for       6  <    |  P 


The  13(8)  function,  which  is  the  response  to  a  Dirac  delta  function 
input  of  a  tracer,  is 


s(e)  =  ^.S(9-^  +  _(1-^ 


o- 


2j3 


^ 


^    38 


=  0 


where  5(6  -  -£)  is  a  Dirac  delta  function  defined  by 


S(8-  |p)=  0 


S  o>o 


for  e  ^  |  A 
for  e  =  |p 


and 


/"s<6-  h)  dO-l 


for  0  >  =  p 


(3-2?) 


for  8  <  |p 


(3-23) 


(3-29) 


2.  Ostwald-de  Waele  Model 

The  Ostwald-de  VJaele  Model  fluid  may  be  characterized  by  the  following 
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relationship  between  the  shear  stress  tensor  and  the  rate  of  deformation 
tensor  (2) 


--H/FI7I)     }s 


(3-30) 


Under  conditions  for  which  a  one-dimensional  rheological  statement  is  valid, 
Equation  (30)  in  rectangular  coordinates  may  be  written  as  (see  Fig.  (1)) 


I   =  -  m 


dV. 


x 


dy 


tf-l  /dV. 


dy 


(3-31) 


dV, 


where  m  and  \)   are  rheological  parameters.  Since  — -  is  always  negative  for 

dy 

the  given  system,  Equation  (31)  reduces  to 


dVv  i/ 


The  momentum-flux  distribution  has  been  found  to  be 


(3-32) 


T   =  fgy  cos  if  (3-7) 

Substitution  of  Equation  (32)  into  Equation  (7)  gives  the  following 
differential  equation  for  the  velocity  distributions 


_  f^  =  (£S2L£°§JL) 

dy       ei 

This  equation  is  integrated  subject  to  the  boundary  condition 


V  -  0 

'x 


at  y  =  S 


to  give 


(3-33) 


(3-3*0 


(3-35) 


where  n  =  -.  The  maximum  velocity  V  occurs  at  y  =  0;  that  is 

in   n  <■  n+1 
m   N   m    '   n+1 


(3-36) 
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Combining  Equations   (35)  and  (36)  yields 

TXW  -  T„  1 1  -   (f  )n+1  3  (3-37) 

The  average  velocity,  V  ,  over  a  cross-section  of  the  film  is  obtained 
as  follows: 


d  0 


-vj'll-cr1]^) 


=  ^v  (3-33) 

n+2  m  \wi 

The  cumulative  age  distribution  of  an  Ostwald-de  Waele  fluid  flowing  through 
an  open  channel  when  the  entrance  effect  is  negligible  is 

F(8)  =  4-  f   V  (y)dy 


y 


=  i  J*  [  1  -  (|)n+1]  d(g)  (3-39) 


'X 


From  Equation  (37)  one  has 

1 


z  =  a  -  V 


1 

n+1 


1    Vx(n+2)  J 
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and 


±  (Z) 


1 


(n+2)  8' 


1  - 


n+1 


(nt2)e 


-ft 

n+I 


(3-^1) 


Substituting  Equations  (40 )  and  (41)  into  Equation  (39)  shows  that 
n+2jn+l  fl3  I    (n+2) 8  J 


(3-42) 


n+2 


Therefore 


f(9)  =  4p 

""     n+2 


0 


n+1   3 
n+2  © 


L    ,n+2v  J 


n  _ 
n+1 


d8 


Wi; 


o 


for  Q  >  SJi 
u  -  n+2 


for  e  <  — 
n+2 


(3-43) 


The  residence-time  distribution,  which  is  the  response  to  a  Dirac  delta 
function  input,  is  obtained  by  taking  the  derivative  of  F(0)  with  respect 
to  0,  that  is, 


E(0)  = 


1 
n+2 


e 


3f 


1  - 


/n+2xnJ 
(nTl)0 


n 

n+1 


for  9  > 


n+1 


—  0 


-  n+2 


f  or  0  <  2±1 
n+2 


(3-44) 


VJhen  the  entrance  effect  is  considered,  the  tracer  is  uniformly  dis- 
tributed in  the  entrance  section.  The  cumulative  age  distribution  at  the  " 
outlet  of  the  system,  which  corresponds  to  the  response  to  a  step  function 
input  of  a  tracer,  is 


Z(0)=|J    dy 


I   J0 
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Therefore 

.0 


F(9)  -  J       1%   ,  .  1 —  dB,  for  e  >  ?il 

— i ] 


.+1  (n+2)?   r         ,  rrr  "  "  n+2 


=  0  for  9<~ 

n+2 

The  response  to  a  Dirac  delta  function  input  is  then 

s(e)  =  JL  . I for  a  >  — 

~     n+2  n  ~  nt2 

^nti;e 

=  0  f  or  fl  ^  — 

nt2 

Comparing  the  distribution  functions  derived  in  this  chapter  for  fluids 
flowing  through  an  open  channel  with  those  derived  by  Fan  and  Hwang  (19)  for 
fluids  flowing  through  a  slit,  it  is  seen  that  the  forms  of  the  velocity  and 
concentration  distributions  are  exactly  identical.  The  differences  arise  from 
the  fact  that  extensive  properties,  such  as  cross-sectional  area,  flow  rate, 
and  total  tracer  injected  in  this  analysis  are  reduced  to  one  half  of  those 
in  Fan  and  Hwang's  analysis.  Thus,  it  is  not  surprising  to  find  that  both 
flow  geometries  result  in  the  same  residence-time  distribution  functions. 
FLOW  OF  THE  BINGHAM  PLASTIC  WITH  SLIP  VELOCITY  AT  THE  T:J3E  WALL 

The  convective  model  for  Bingham  plastics  will  be  treated  by  using  a 
generalized  velocity  profile  for  which  the  slippage  at  the  tube  wall  is  con- 
sidered. The  momentum  flux  distribution  for  the  flow  of  fluid  through  a 
circular  tube  has  been   given  as  (2)  (see  Fig.  (2)) 

Tra-(-?0~2LpL>r  on 
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in  which  P  =  p  -  pgz  and  where  p  is  the  pressure  acting  on  the  given  system. 

The  one-dimensional  rheological  statement  of  the  Bingham  plastic  in 
cylindrical  coordinates  is  found  to  be   (2) 

dV„ 
T       =  t     -  u„  — £;  for    "C       >   T  (3-Li-6) 

crx         0      r0  dr  **  -       0  w       ' 

Combining  Equations  (^5)  and  (^6)  and  applying  the  boundary  condition 

V.=  Vp  at  r  —  R,  the  velocity  distribution  is 
x    it 


d-io) 


where 

°"~  % 

and  V  is  the  maximum  velocity  of  the  fluid.  The  average  velocity  of  the 

fluid  is  calculated  by  summing  all  the  velocities  over  a  cross-section  and 

then  dividing  by  the  cross-sectional  area  as 

2   R 
J   /  Vx(r)rdrd9 

J   J   rdrdB 
0   0 

0  (i  -l0r 
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=  y«  +  (1  -  \)  vR  (3-50) 

■where 

ol  =  3  o   ?  o  =  i(3  4  2  ^  +  ^j  (>5l) 

The  minimum  residence  time  of  the  fluid  will  be  defined  as 


°  min   v  ~       v 

m         m 


o( 


(3-52) 


2(l-(l-|)v] 

where  V  is  a  slip  velocity  factor  and  is  defined  by  V  =  vr/Vx*  ^  exPlici't 

expression  of  £  in  terms  of  V  can  be  obtained  from  Equation  (48)  as 
R  X 


|  =  L  +  (1 :- $  )  /  5l1-!s  (3-53) 

m   VR 

Substituting  Equation  (52),  V  =  vr/^x  and  Vx/Vx  =  0  into  ^uat10*1  (53)  yields 


r   t        ,      2  1-  l-^yl-cj/e 

I  =  L  +  (1  -  1 0)  J  -* ^— I (3-54) 

R    °        °  N        2(1  -V) 


and 


j 


&(f)-   -  __j^  (3-55) 

46 


Jl-/J[l-  (1-|)4<A/2S 

In  response  to  a  unit  step  function  input,  the  cumulative  age 
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distribution  is 


?(0)=  -4-  I       Vx(r)rdrd0 


(3-56) 


-  i^w  +  z2(0) 


(3-57) 
(3-58) 


where 


h.™  =  ir  1 


*V 


.-,   0 


V  rdr 
n 


!^(i-a.|)yj.Us{e-- 


l  -  (l  .  |)V] 


l 


and 


z2ce> 


x   0 


fN^-vl1-?^!] 


rdr 


u-U' 


(3-59) 


^/^^-'"^■^KJ^lU   (3.60) 


Substitution  of  Equations  (5*0  and  (55)   into  Equation  (60)  reveals  that 


Z2(8)=J 


G 


0 


d> 


2[l-  (1-|)V) 


i  .  (i  -<*)V-  A 

2     29 


.+ 


i-io  i  a-$0)°« 


a/T-V 


0-1  d6 


JT^v 
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./iv  I      .rrz^z      "-'■",.  /T7«>    '  3*        a 3/2 


/T^7  I      VI~^2 


^ii^2      3*         e 


^[i-u-  %)i%    A/2[i.(i.  |)yje-d} 


3^ 


0 


1/2 


for 


a 


_  ^  e  <  i     (3-61) 

2[i-(i-  f)v]        r 


The  complete  expression  of  F(9)  is  then 


F(9)«!kiuj. 


* ,  ,  (i-$o>°< 


f  (l4n)A2         1-J( 


2V      ,71^7 


//i^V<*      Wi-Ve2 


1 — —    rrr~  + 


34 


0 


3/2 


3eC 


,1/2 


=    0 


^    &<&£$ 


*°*     **& 


(3-62) 


=  1 


for    ©  >  i 


where 


A=l-  (i-|)V 


(3-63) 


The  residence-time  distribution,  which  is  the  response  to  a  Dirac  delta 


function  input,  is 


E(e) 


2V^<,*      ol.^^io^    f      ^0  ^ 


_  "■*<) 


s<e-£>+ 


(' 


) 


27 


for  —  ^9  ^  i 

2A     ^  (3-6*0 

=  0  elsewhere 

When  the  entrance  effect  cannot  be  neglected,  the  F- function  can  be 
calculated  by 

-i   r2TT  rr 
F(0)=  -i-  I       rdrd0  (3-65) 

Thus 

^  S         2A         iJlIV     l  <*  29A/I37       oJl~7J 


f0P    f^0"      J 

-  °  for    e  <  ^ 

=  1  for     e  >  ~ 

y 

and  the  corresponding  S- curve  is 

°^  2A  / _2  I  / r       n~z  J 


(3-66) 


/it/ 


(3-67) 


=  0  elsewhere 

When  y=  0,  the  slip  velocity  at  the  wall  of  the  tube  is  zero,  Equations 
(62),  (6^0,  (66)   and  (67)  respectively  become 


a/^1q     J  29-o(      ^J2  *Q     ^29-01  I 

301    63/2       3^2      e1/2  J 
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for  0  >  - 
~  2 


=  0  for   6^ 


2  , 

5(6>  =  %  S(6-  f )  +  i^  .(-k=  +  l.|0) 


=  0  for  0  < 

2 


m ,  %!  Us(9.  |) + (1.Jo),|fkS.  hk + Lk  ] 


=  0  for  e   <  | 


2  L 

so)  =  %  $®.  %)  +  id^d^ + i-i) 

*       2      262   ^     0 

for  9  >  ^ 
-  2 

*  0  for   0  <  - 

2 


(3-68) 


for  e  >   |  (3-69) 


for   0  >  |  (3-70) 


(3-71) 


When  V=0  and  jjQ  —  0,  the  above  distributions  reduce  to  Equations  (2-1),  (2-2), 
(2-3),  and  (2-4)  respectively,  -which  have  been  derived  for  the  Newtonian 
flow  (4,5). 

A  family  of  numerically  computed  F  and  E  curves  with  flow  behavior 
index  £Q  and  slip  velocity  factor  V  as  parameters  are  shown  in  Figures  (3) 
through  (12).  Since  molecular  diffusion  has  been  neglected,  certain  points 
of  discontinuity  may  occur  in  the  distribution  curves  due  to  the  irregular 
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Theological  behavior  of  a  Bingham  plastic  fluid  and  the  slip  velocity  at  the 
tube  wall,  l&en  £,  becomes  0,  the  flow  system  reduces  to  the  Newtonian  case 
and  the  resulting  concentration  distribution  are  analogous  with  those  proposed 
by  other  investigators  (4,5).  However,  the  E-curve  and  F-curve  are  deformed 
in  actual  observations  since  the  effect  of  molecular  diffusion  does  exist  in 
actual  systems. 
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IV.  DISPERSION  MODEL 

A  generalized  mathematical  expression  for  diffusion  in  the  flow  system 
at  constant  temperature  and  pressure  has  been  given  (2)  where  D,  the  molecular 
diffusion  coefficient,  is  assumed  to  be  constant.  Suppose  that  there  is  no 
chemical  reaction  occurring  in  a  fluid  which  flows  through  a  circular  tube 
and  that  the  concentration  distribution  is  symmetrical  about  the  central  line 
of  the  tube.  The  diffusion  equation  becomes'  (see  Fig.  (2)) 

£C  _  _  ,&2C  ,  ldC,  S2C,   „  ,  .  £C 

it  -D^2  +  r  S7+  ^2>  -  Vx^)  a3E  (4-1) 

Taylor  (8)  has  shown  that  under  the  condition 

E  >^  D  (4-2) 

the  transfer  of  C  along  the  tube  by  molecular  diffusion  is  small  compared 

&    C 

with  that  produced  by  convection  and  thus  the  — — ^  term  may  be  neglected.  He 

has  also  shown  that,  when  the  time  necessary  for  a  radial  variation  in  C  to 
die  down  owing  to  radial  diffusion  is  much  shorter  than  the  time  necessary  for 
an  appreciable  change  in  C  to  occur  through  longitudinal  convection,  an  approxi- 
mate solution  can  be  obtained  after  certain  simplifications.  The  model  corre- 
sponding to  such  a  solution  of  the  partial  differential  equation  of  diffusion 
is  called  the  dispersion  model. 

Taylor's  assumption  can  also  be  expressed  as  (6) 

Vm    (3.8)2D 
This  is  applicable  to  all  types  of  flow  models  as  long  as  fluids  flow  through 

circular  tubes  and  the  ~  term  can  be  taken  as  zero  in  the  derivation  (6). 

dX 

If  rectangular  coordinates  are  used  (see  Fig.  (13) )»  a  corresponding 
simplified  expression  of  the  partial  differential  equation  for  diffusion  is 
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In  order  to  find  the  conditions  under  which  Taylor's  limiting  condition  may- 
be valid,  it  is  necessary  to  calculate  how  rapidly  concentration  becomes 
uniform  due  to  molecular  diffusion.  The  governing  differential  equation  is 
2 

du    *y2 

subject  to  the  boundary  conditions 

—  =  0  at     y  =  0  and     ±  %  (4-6) 

sy 

It  should  be  noted  that  the  velocity  distribution  of  the  fluid  over  the 
cross  section  does  exist.  Introducing  the  diraensionless  variables 

tV    j 

one  obtains  Equation  (5)  in  diraensionless  form  as 

(4-7) 


3C  _  Dt      a2C 


The  technique  of  separation  of  variables  is  used  to  solve  this  equation,  or  in 
other  words,  C  =  F(9)G(4).  Substitution  of  this  relation  into  Equation  (7) 
gives 

Fde     s2G  de2 

This  equation  has  the  following  solutions: 

F(6)  =  Ax  e"a  G  (4-9) 


i  2  ?  I  2  2~ 

G(0  =  A2  cos  (J^i)  +  A3  sin  (J  ^4- 4)  (^-10) 
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-I»  ^2  an<^  ^3  are  constants  ^°   Dc  determined  by  the  boundary 
conditions.  JU  must  be  zero  because  the  concentration  distribution  is 

.etrical  around  the  center  line  of  the  tube.  Equation  (6)  requires  that 


i  2  ° 
sin   /  £JC   =  0 

a/  Dt 


Therefore , 


/^X  =  r.TT.     n  =  0,  1,  ....  n  (4-11) 

n  Dt 

The  root  of  this  equation  corresponding  to  the  lowest  value  of  d,   is 

<C  =  n2  4  (4-12) 

1  r 

and  the  solution  of  Equation  (7)  takes  the  form 


C  =  A  e~V  cos  (  JSLL.  £)  (4-13) 

w  lit 

The  fact  that  the  concentration  at  6  =  0  and  £  =  0  is  1  ensures  that  A  =•  1. 
From  Equation  (13 )•  the  time  necessary  to  decay  the  central  concentration 
down  to  e~  of  its  initial  value  is 


(4-14) 


Jr 

Dlf 

=1  u, 

101 

D 

i  implies 

that 

L 
Vm 

>>  0. 

101 

D 

(4-15) 

in  order  that  Taylor's  second  limiting  condition  may  be  applicable.  Equation 
(15)  means  that  the  time  necessary  for  appreciable  effects  to  appear,  owing 
to  convective  transport,  is  long  compared  with  cross-sectional  variations  of 
concentration  through  the  action  of  the  molecular  diffusion. 

The  systems  studied  in  this  chapter  are  those  in  which  the  dispersion 
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in  steady  flow  is  due  to  the  combined  action  of  convection  parallel  to  the 
plates  and  molecular  diffusion  in  the  direction  perpendicular  to  the  plates. 
Restrictions  suggested  by  Taylor  are  basically  followed  in  the  present  analysis, 
FLOW  THROUGH  A  SLIT 
1.  Bingham  plastic  model 

The  diffusion  equation  becomes 

■I!  »  D  i-C   v(y)  -^  (4-4) 

Vfoen   a  one-dimensional  rheological  statement  is  valid,  the  velocity  distri- 
bution of  the  Bingham  plastic  fluid  in  rectangular  coordinates  is 


(1       >  ^ 


(-y1=\[1-^ ^  ].  for  Irl  S|y0|  (3-14) 


(1  -  i0) 


where 

2 


V  =  PZ*     cos^  (1  -  I    )' 


^0 


(3-13) 


and 


vx=  3  f\  (3-1?) 


Combining  Equations  (4)  and  (3-14)  yields 

2.        ,       ttZ\     ;  \2 


g-^-^j^ji,      for|yl,|7ol 


^  =  D  a!s. v 


(4-18) 


iT  =  u  73"v»di  fOT  iri  <|yfl| 


Introducing  the  dimensiohless  variables 


k5 


tv 

ft  —    x  _  t         >     I  V  I      .-   X 
T.  b      I  5.        I    L 


t  I  % 

one  obtains  Equation  (18)  in  dimensionless  form  as 

2 


ac  =  Dt   a  C   J3_ [■]_ 

^   S2  '  zf"  2M 


(t-4nri 


0' 


(l-^o)2j 


_  Dt  a>  c   3  a  c 
£2'  a^2"'  2?'  M 


with  the  boundary  conditions 


-£■£     forj   I 


for  ^  <  £ 


H 


=  0 


at  £  =  0  and  1 


(4-19) 


(4-20) 


since  the  plates  are  impermeable.  For  convenience  the  derivative  following 
the  mean  speed  of  flow  is  introduced.  Relative  to  these  axes,  the  velocity 
distribution  is 


V  ,  =  V    -  V 
xl       vx        x 


-l$-« 


1  (^-*o>2 
"2p  a-o2 


] 


V. 


=     %-^\ 


*»    $  >l 


£°*  4<  40 


(4-21) 


or 


Vxl= 


V 


x 


(*0+    2)    (l-40)' 


((l^o)3-3(^^o)2J     for   J    >l 


(4-22) 


a 


&-i0>- 


*°n  <i0 


'0  •    "' >0' 

Let    7i=  7-  9    ■  Equation  (19)  becomes 
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E£  -  £l  Jl£  _      i f  n-M3  -  3^-  n2 1  * c     for  >  > ; 


2LV      70'  "^    '0'    J    3fl,  ^  -  s0 


(^-23) 


«  St  ^fc 1 (1  „  >    )3  _^C_  for    <  <  £ 

Since  the  moving  axes  which  move  with  mean  fluid  velocity  are  introduced,  the 
transfer  of  tracer  concentration  across  the  plane  at  which  ^  is  constant 
depends  only  on  the  cross-sectional  variation  of  tracer  concentration.  A 
partial  equilibrium  may  be  assumed  for  a  small  cross-sectional  variation  in  C, 
If  _iL2  is  considered  to  be  independent  of  £',  the  small  cross-sectional 
variation  in  C  can  be  calculated  from 

*2  H2   (^+2)(l-0      °         °    »7l  ° 


(i0+2)(l-^)2      "   **l 

The  concentration  distribution  is  found  as 


(I-^O)3-^-        for    4  <  i 


(4-24) 


2(4Q+  2)(l-^0)2^      ^1 


for    £,  >  £ 


0 


(4-25) 


.2 

=r  cn+ 1 ~K  (-££-)  (i  -£J  42     for  >  <<> 

0       2(40+   2)(l-t0)2Dt      ^i  *°'    *  *        S0 

The  rate  to  which  C  is  transported  across  a  section  at  */-,  is 

Q  =  2S://   CV  ({)  d£  (4-26) 

Therefore 


17 


^=I°(c0+M(i-i0)V]«H)3d* 


+  /  [co  +  M  f  CH>>Y"  |«-«0^J  N((i-i0)3-3(4-^o>2j  d<      C*-2?) 


v.-here 

,2 


Ha 


Dt      ^Y1     2(^+    2)(l-*0)' 

^^ 

(£0  +  2)(i-40)2 


and 

JL  =  -  -L  (i  +  22  z   +  a  >  2}  (H)2  <^x_  J^ 

2SW    105     16  0   16  *0'  v£  +2    DL    £  >?, 

0  -L 


(4-23) 


Since  the  cross-sectional  variations  in  C  are  assumed  to  be  small  compared 
with  those  in  the  longitudinal  direction,  ^CwiA7i  ^s  near-*-Y  ^e  same   as  ^Cm/^l 
if  Cjjj  is  the  dimensionless  mean  concentration  over  a  cross  section.  Equation 
(23)  may  be  written  as 

Q  =  -2$wJL  (1+  33  /  .  21  >2n  ^""V2  $>\  *Cm  (^29) 

Comparing  this  with  Fick's  lav:  of  diffusion,  it  can  be  seen  that  C_  is 
dispersed  relative  to  a  plane  which  moves  with  velocity  V  with  a  diffusion 
coefficient 

s  =  -I-  (i  ♦  22  /  +  ZLf  )  (ilk)2  Li  (MO) 

105     16  0   16  0    £Q+2    D 
which  is  called  the  effective  or  apparent  dispersion  coefficient. 

It  should  be  remembered  that  two  conditions,  Equations  (2)  and  (15).  are 
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necessary  for  the  given  approximate  solution  of  the  equation  for  diffusion  to 
be  valid  to  interpret  the  longitudinal  dispersion  of  a  solute  in  a  stream. 
From  the  first  condition  and  Equation  (30),  one  knox-rs  that 

105  (   16  *0 +  16  V  V1  ~  >   D  (4"31) 

in  order  that  the  longitudinal  molecular  diffusion  may  be  negligible  compared 
with  the  dispersive  effect  represented  by  S.  The  second  condition  ensures 
that  the  time  necessary  for  a  cross-sectional  variation  in  C  to  die  down  owing 
to  cross-sectional  diffusion  is  much  shorter  than  the  time  necessary  for  any 
appreciable  change  in  C  to  occur  through  longitudinal  convection.  From  this 
condition  and  the  fact  that 

V  =  -r-2-.  VY  (4-32) 

m   £0+  2  x 

one  obtains 

JL  >>  0.034  (£0+  2)  JL  (4-33) 

vx 

Combining  Equations  (31)  and  (33)  gives  the  conditions  under  which  the  given 
solution  is  valid. 
2.  Ellis  Model 

It  has  been  shown  that  the  momentum  flux  is  related  to  the  velocity 
gradient  according  to  (2) 

X     =  (  ° "  ?L)y  (4-3*0 

yx     l 

and  the  fluid  which  has  the  shear-stress  and  shear-rate  relation 

-fs.-^  +  <f_|-r  I11"1)  X  (4-35) 

dy     J  0    J 1 1  yx '      yx 
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is  called  the  Ellis  model.  y>   ^>  and  n  are  measurable  parameters.  The 
Ellis  model  is  the  simplest  and  most  generally  useful  three-constant  model. 
It  describes  properly  the  lower  limiting  viscosity  >J   .  This  model  seems  to 
have  sufficient  flexibility  to  fit  the  data  for  various  types  of  fluids. 
Substituting  Equation  (3^)  into  Equation  (35).  one  obtains 
dV. 


-iM^it^r1)^ 


=y0  C=£*)y  +  ^  (^)n-  y*  (4-36) 

This  equation  can  be  integrated  subject  to  the  boundary  condition 

Vx  *  0  ,  at  y  =  1  S 

The  resulting  velocity  distribution  is 

\  '   ^0  Pf)  4  (l  -  (f)2j  +  *,  (=f  >*  -^  ll  -  (f)n+1](^-37) 


The  maximum  velocity  is  the  velocity  at  y  =  0,  that  is, 


V    =  if0  (=££)  i!  +  tf     (rAP)n  _*^ 
B  °        I        2  1       L  ntl 


(4-38) 


=  \l  +  V 

Expressing  Vx  and  Vx  in  terms  of  V^  and  Vn£»  one  obtains 

V  =  %  V  -  +  £ii  V  (4-39 ) 

and 

Combining  Equations  (40 )  and  (4)  yields  the  governing  partial  differential 
equation,  that  is 
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"  *2  h2   v  »*      ?       m2    s     J  (ar 


58 


(4-41) 


with  the  boundary  conditions 


ac 


=  o  , 


at   £  =  0  and  1 


(4-42) 


If  one  defines  concentration  and  velocity  relative  to  axes  which  move  with 
the  mean  velocity  of  fluid  flow,  the  velocity  distribution  is 


Vxl=Vx-Vx 


-V§-r>+  V*2(nk-*    > 


(4-43) 


and  the  transport  equation  becomes 


2  2  n+1  -, 


The  distribution  of  C  under  the  conditions,    <S>C/©>8  =   0  and    ^c/^'/-i  =  constant, 
is  found  to  be  in  the  form 


c = c°+  J-  M  ft v^3  •  s  v^+  v*2  ti& "  .^£3). 

■*■   '1 


■]](«5) 


where  C  is  the  concentration  of  the  tube  at  2,  =  0. 

The  rate  at  which  C  is  transported  across  a  section  is 


=  &!$J 


CV^d* 


(4-46) 


Inserting  values  of  C  and  V  -,  into  Equation  (46),  Q  is  found  to  be 
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Q=     ^  ("^H^  ^+  27<n+4>(2n+5>  ^ 


l6(n+l?  (n+?) Vnl  Va2 1  (4-47) 

405(n+2)(n+4)(n+6)       J 

Since  BC/^yJ1   has  been  assumed  to  be  independent  of  £  ,  C  may  be  replaced 
by  C  •  The  rate  of  transfer  of  matter  in  a  slit  due  to  a  diffusivity  2  is 


2^3  f^a  (4JW) 

L    ^1 


Comparing  Equations  (47)  and  (48),  it  can  be  seen  that 


2 '  °  v^+  stC^W  v^  +  w(l§>B&  ^  v]  (*-*?) 


2  =  11.8  v2_  j. 


945 

It  should  be  recalled  that 

S  >  >  D 

in  order  that  the  axial  molecular  diffusion  be  negligible  and  that  the 
following  condition  must  be  satisfied  so  that  Taylor's  second  assumption 
may  be  true: 

■L  s.  >  o.ioi  -i. 


where 


When  if,  is  zero  the  Ellis  model  reduces  to  the  Newtonian  case  and  the 
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effective  dispersion  coefficient  becomes 

because  Y  0  becomes  zero.  For  the  case  4>  =  0,  V  ,  becomes  zero,  that  is, 
m2  J  0      ml 

for  an  Ostwald-de  Waele  model,  the  effective  dispersion  coefficient  is 

E=3(n-|4)(2n+5)  "D^  <^2) 

FLOW  THROUGH  CYLINDRICAL  TUBES 

If  the  transfer  of  a  tracer  along  the  tube  by  the  molecular  diffusion  is 
small  compared  with  that  produced  by  convection,  the  transport  equation 
becomes  (see  Fig.  (2)) 

The  boundary  conditions  are 

—  =  0  at  r  =  0  and  R  (^-53) 

because  the  concentration  distribution  is  symmetrical  around  the  tube  and  the 
wall  of  the  tube  is  impermeable.  Introducing  the  following  dimensionless 
variables 

t   x     u    r         x 


Equations  (1)  and  (53)  become 


5C   Dt  /  a  C  ,  1  ±Cs   _  \  ^C  (i|-5^) 


*e    r2  w   s  n    vv  ^7 


'X 


53 


ana 

JL2  =  0  ,         at  NOandl  &-55) 

The  condition  under  which  Taylor's  approach  may  be  valid  can  be  expressed 
as  (6) 

Jl  >>  _2_  (4-56) 

Vm      (3.8)2D 

Velocity  distributions  which  can  be  represented  by  the  Bingham  plastic 
and  Ellis  model  will  be  considered  in  this  section. 
1.  Bingham  Plastic  Model 

The  velocity  distribution  of  Bingham  plastic  fluid  has  been  derived 
as  (see  Equation  (3-^8)) 

as  V  for  fc  <:  L 

m  *   ->0 

and 


(^-57) 


;;re  are  now  considering  the  convection  across  a  plane  which  moves  with  the 
mean  speed  of  flow.  Writing 

fi-T-fl  (^-59) 

the  velocity  distribution  relative  to  these  axes  is 


V  ,  =  v  -  V 
xl    x    x 
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($ .  i)  -  s  £kio>!] 


<* 


d 


v. 


0-  -  u 


X 


and  Equation  (54)  becomes 


for    |    >     ) 


for     ^    <    1 


(4-60) 


2 

a_C  _  Dt   /  ^   C       1  ^Cn 


d^2   i  **    l<* 


(4-D 


3^    forU>0 


Dt 
R2 


a2c 


1    £C< 


C — 2  +  7  TT' 


(2  .  d  ^c 


for  *  < 


(4-61) 


The  small  radial  variation  in  C  can  be  calculated  from  Equation  (6l)  by 
letting   ac/<3>0  —  0  and  considering    ac/^>/-i  ^°  ^e  independent  of  £,  that  is, 


Dt   /  £   C       1    SC) 

r2  *i2  +  s  n 


=  Dt  (1  _±_  (t±C\] 
r2  L'g  »|    v3  *fc'j 


I 

(H0): 


jL2 
^(1-V2-1  ^1 


=    (2-i)  -i-2 


for 


lor 


i*l 


t<t 


(4-62) 


Equation  (62)  has  the  boundary  condition 


ac 


=  0 


at     I  =  0  and  1 


The  solution  of  Equation  (62)  can  be  obtained  as 
2  rS0  ,    f*0     , 

0  =  0+^(^2.)/      Id   J      (i-Djdj 

0       Dt     ^7i       JQ       I      >  JQ      &  * 


(4-63) 
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R-    /^C 


Sid*  r^-i)-2<iiklf 


Dt    d7i  4  $<sl  to)  o>  (i-|0)2Js   5 

=  cn+£  (±2-)(l  idJ    (i-l)U  I 


Inserting  c*  =  =  (3+2J0+|Q)  into  Equation  (6k)   yields 


for  £  ^  ^0 


for  $  <  £ 


(4-64) 


C  =  C0  + 


ac 


Dta(l-I0)' 


s0    ^  9  *0  ^  8  $ 


2    ^7-l  L  4       3  ^0       12 


+  12t    +  ir  ln-il 


*>*  *  >  t 


=  c0  + 


?,  _£C 

DtcJ(l40)2   ^7 


4  1,2 


^c  ri _  2 1     i^  _  i_r ir 

^U       350       2^0         123oJ^     ,   for   |  < 


(^-65) 


where  CQ  is  the  value  of  C  at    \  =  0. 

The  rate  of  transfer  of  C  across  the  section  at  ^  ^s 


Q  =  2ITR 


i      C\lld^C  CVxl  <J)*d! 


(4-66) 


'0   ' 


Substituting  Equations   (60)  and  (65)  into  Equation  (66),  the  rate  of  transfer 
of  C  across  the  section  at  ^     is  found  to  be 


Q  =  - 


4-9 


'2  .,  44 


+ 


l6.2      ,4 


2L3(3t2^)2(1^0)U8       35     °       15     ° 


$0     15  5o      5^0 
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The  rate  of  transfer  of  matter  across  a  cross  section  with  the  diffusivity 
E  is 

q  *  .  EO  (±£B)  (4.68) 

L    ^3/ 

It  should  be  noted  that  <^c  /51/   is  indistinguishable  from  «5c/^7-j  as  long 
as  ^C/^7  ]_  is  independent  of  \.     Thus  it  will  not  be  unreasonable  to  define 
the  effective  dispersion  coefficient  E  of  Bingham  plastic  fluids  flowing 
through  a  cylindrical  tube  as 

_2-2 
E  _  R  vx  f  3      t£  fc   +  I6fr2+ V4      28c5      2v6 

*  2D(3+2|o+^)2(l4o)U8  "  "35  °   15  °  l°"15     °  "  5  S° 


For  the  case  jn  =  0,  Equation  (69)  becomes 

E  =  -^SF  (^70) 

This  agrees  with  the  result  for  the  Newtonian  flox*  derived  by  Taylor  (6). 
Since  the  final  expression  for  S  appears  to  be  peculiar,  it  may  be  desirable 

to  evaluate  the  value  of  E  at  5  =  1.  As  $0 «-  1,  In  ^n  may  be  expressed 

as 

m  ^0=  -  (i40)  -  I  (i-lo)2  -  \  a-^o)3  -  I  (i-!o)4-...    V*r7i) 

Combining  Equations  (71)  and  (69).  it  can  be  seen  that  the  effective  dispersion 
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coefficient  S  becomes  zero  as  |0  — «-  1.  This  is  reasonable  since  we  have 
neglected  longitudinal  diffusion  in  the  analysis. 

The  conditions  under  which  the  given  approximate  solution  can  be  used 
to  interpret  longitudinal  dispersion  of  a  solute  are,  as  stated  previously* 


£ 


(4-72) 


where  S  is  given  in  Equation  (69)  and 

(3  +  2|  +  &L  R2 

Vn       6vx  (3.8)2D 


2.  Ellis  Model 

The  momentum  flux  distribution  for  flow  of  a  fluid  through  a  circular 
tube  is 

P  -  P 

Tw  =  (— -)  r  (4-7*0 

rx   v  2L 

where  (P^-Pr  )  ^s  "^e  result  of  a  pressure  gradient  and/or  gravitational 
acceleration. 

For  the  case  when  a  one-dimensional  rheological  statement  is  valid, 
the  relationship  between  the  shear-stress  and  shear-rate  written  in 
cylindrical  coordinates  is  (see  Fig.  (2)) 

*VV   ,,    .  ,    ,n-l 


x  _ 


dr 


=  <*o  +  *iK*l    >t«  <*-«> 


where  •j*   ;  ^ ,  and  n  are  parameters  to  be  determined  by  experiment, 
Substitution  of  Equation  (74)  into  Equation  (75)   gives 


dVv   ..   ,.  Aov    .  ,_ap  *n 


-  ^  =  J  (=*£),  +  y  (-"^r)  (4-76) 

ar     0  2i»     j1  2i_. 
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The  boundary  condition  that  V    =  0  at  r-R  can  be  used  to  evaluate  the 
integration  constant.     Then  the  velocity  distribution  becomes 

The  maximum  velocity  occurs  at  r  =•  0;  thus  it  has  the  value 

-AP<4>      _2           APn^Rn+1 
v  =  ( 2)  £-  +  (=££)    li (4-78) 

vm       *>   2L     '    2    T  v  21/         n+1  v      r   ; 

=  \l+Vm2  ^-79) 

Combining  Equations  (77)  and  (79)  yields 

v-*-Ji-<i>2J+%2[i-®n+1]-  (*-«» 

The  average  velocity  Vv  is  calculated  by  summing  up  all  the  velocity  over 

Jib 

a  cross  section  and  then  dividing  by  the  cross-sectional  area,  -which  yields 

V  =  1  V,  +  S&  V  (4-81) 

x      2    fill       nt3    m2 

The  mathematical  expression  of  the  diffusion  equation  for  an  incom- 
pressible fluid  flowing  in  a  circular  tube  with  the  assumption  that  the 
axial  diffusion  may  be  neglected  is 

4f=M^+^)-Vx||  (,-83) 

It  is  convenient  in  this  discussion  to  define  concentration  and  velocity 
relative  to  axes  moving  with  the  mean  flow.  Relative  to  these  axes,  the 
velocity  V  -,  is 
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V  -  V  -  V  -  V ■  f  1  -  (|)2]  +  V   I  2  .  (|)n+;Ll         (4-83) 

xl    x    x    mll«c    R  J    n2^  n+-3    R    J 


Substituting  Equation  (83)  into  the  governing  partial  differential  equation 
and  writing  it  in  dimensionless  variables,  one  obtains 


n+lN  1  £  C 


C  _  Dt  ,  a  C  .li.Cv   1  (,r  A   «.  2v  ,  ..  ,  2    .  n+lv  1  _£C 


(4-84) 


The  distribution  of  C  in  the  case  when  &  C/«3  &  =  0  and  -— ~-  is  independent 

d"l 

of  \   is  found  as 


c  -  c0+  -fL  (|f-)(Vnl<|  f-  £  k*>  +  ^(|  f -  £  i*>)l 


(^85) 


by  applying  the  boundary  conditions  that  £  C/^  £  =  0  at  ^  =1. 

The  rate  of  transfer  at  which  C  is  transported  across  a  section  at  7n  is 


9    r1 


(4-86) 


Replacing  C  and  V  .by  Equations  (85)  and  (83),  it  is  found  that 

q=2.tr2     K2    f^c,   f\v     ,is2      1    iV\2,n2       1   tn+3  I 


=  -  2TTR 


2       R^     ,  ^C 


C4*-> 


(n+1)' 


Dt  Vx     *Ti  1384  4(n+3)3(n+5) 


V. 


ra2 


(Yl+l)(ri4ll)  "J  /a  ooN 
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Comparing  this  with 

Q  =  -ttr2  1    ml  (4-88) 

L      d71 

the  effective  dispersion  coefficient  E  is 

D  1 192  V^+  2(nt3)3(nt5)      12(n+3)(n+5)(n+7)  ^  m2J  ^  ^ 
It  should  be  noted  that  the  conditions 

E  =>  >  D  0-90) 

and 

L     >;>    R2 


\l+Via2      (3.8)2D  ^91) 

must  be  satisfied  in  order  for  the  analysis  to  be  valid. 

When  ^  n  =  0,  V  _  becomes  zero  and  the  Ellis  model  reduces  to  the 
Jl      m2 

Newtonian  case  and  the  effective  dispersion  coefficient  becomes 

2  -2 
R  V 

E  -  tst  e*-92> 

If  y   =■   0,  that  is,  for  the  Ostwald-de  Waele  model,  V-  is  zero  and  the 
effective  dispersion  coefficient  becomes 

R   V 

E  "  2(n.3)(n^5)  IT  (^93) 

which  is  identical  to  that  derived  by  Fan  and  Hwang  (19)  for  the  Ostwald-de 
Waele  fluid.  The  limiting  condition  for  the  Ostwald-de  Waele  model,  as 
derived  by  Fan  and  Hwang,  are 


2  -2 
R  V7 


2(n+3)(n<-5)    D 
and 


2  >  >  D  (4*94) 
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Vx   (ntl)(3-8)TD 

The  equation  of  continuity  for  C  in  the  given  process  is 

^Q  ,  =  _  ay  — -2  (4-96) 

in  which  A  is  the  cross-sectional  area  of  the  system.  Substituting  for  Q 
from  Equation  (38)  the  dispersion  relative  to  axes  moving  with  speed  Vx  is 

governing  by  the  equation 

a  C     .   a  2CL 

Since  the  concentration  distribution  has  been  reduced  to  one  dimensional 
statement,  one  can  replace  C  by  C  in  Equation  (97) »  that  is, 

This  has  to  be  solved  subject  to  the  boundary  conditions 

B.C.I.     C=0     for   8=0     and    Y[  *j*  0         (4-99) 
B.C.  2.     C  =  6(8)  for   6=0    and    *J  =  0         (4-100) 
B.C.  3.    f°°C  dl  =  1  (4-101) 

The  solution  subject  to  the  given  boundary  conditions  is  (20,  21) 


C  = 


1  *? 


c 


2 


Ti(H~)e 


'XJ 


172 


exp 


2 


^  VYL  J  VVL 


::A 


Since  ^  =  f  -  ft  has  been  defined  in  this  analysis,  one  may  replace  1-,  by 
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*[  -  0,  which  gives 

C  =  -p 1 375  exp  f-  -Cffi)..  ]  (4^103) 

1  VXL  J  VXL 

An  experiment  designed  by  setting  »7  *  1  yields  the  outlet  tracer  concentration 
distribution  function  to  be 

C=  E(9) 


exp 


(l-9>2 


1  V-  J  V- 


(4-104) 


where  E/VJL  is  a  function  which  varies  with  dimensionless  parameter  ft  vZ/D 

2—2 

(or  R  V^D)  and  the  flow  behavior  index  describing  different  flow  models. 

Integrating  the  given  distribution  function  from  0  to  9,  we  can  obtain  the 
concentration  distribution  corresponding  to  a  step  function  input,  that  is, 

C«  F(B) 


9  ? 

:  f  1 _  exp  (  -  Sk&L  1  d0  (4-105) 

J°  zk(J-)6)l/Z         l     4(^)6 J  . 

1   V_T.   J  V  T. 


A  family  of  E(0)  and  F(0)  curves,  which  have  been  numerically  computed  and 

shown  in  Figures  (14)  and  (15) »  are  plotted  with  E/V^L  as  parameters.  The 

-  2-2       2-2 

correlations  among  E/V  L,  the  dimensionless  parameter  $  Vx/D  (or  RyD)  and 

different  flow  behavior  indexes  are  shown  in  Figures  (16)  and  (17). 

When  E/VL  is  very  small,  the  value  of  E(0)  are  essentially  zero  except 
for  the  values  of  9  close  to  1.0.  For  this  small  E/VXL,  the  residence  time 
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0 
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Di.TSGnsionlecs      TirrsG  ;   9 
Fig.     14.     Gsneralizsd      E- curves     with 
G3    the    parameter. 
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distribution  functions  can  approximately  be  reduced  to 


3(9)  = 


Jirfcrr)] 


8XP 


■which  is  a  nornal  distribution  function  with  mean 
and  variance 


cS  -  2C-&-) 

VXL 


(4-io6) 


(4-107) 


(4-108) 


The  approximate  cumulative  age  distribution  can  therefore  be  expressed  by  (6) 
F(8)=  |  -  1  erf  [  |  (l-0)(E/VxL)*]       for  6  <  1 


(4-109) 


=  l  +  lerf[|(9-l)(E/VxL)lJ 


for  6  >  1 


where 


erf  z  =  2TT 


*L 


z    ,2 

e"Z  dz 


If#  one  differentiates  Equation  (104)  with  respect  to  8  and  sets  the 
derivative  equal  to  zero,  the  solution  shows  that  the  maximum  value  of 
E(8)  occurs  at  time  Q  ,  that  is, 


m' 


or 


8„  +  2(=E-)9  -i=o 


at    8=8 


m 


(4-110) 
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Therefore 


,  m 


t 


E  v2 


L  V 


*   /B 


(—-)  +  1  2  -  (—-)  (4-111) 


V- 


or 


((JL)2  +  il1/2+  (JL) 

This  equation  shows  that  the  maximum  point  is  close  to  1  when  E/V.JL  is  far 
smaller  than  one.  The  value  of  the  residence  time  distribution  at  the 
maximum  point  is  obtained  by  substituting  Equation  (111)  into  Equation  (104). 

The  stimulus-response  analysis  done  by  Taylor  provides  a  convenient 
mathematical  technique  for  finding  the  age  distribution  of  a  tracer  passing 
through  closed  vessels.  The  usual  method  of  finding  the  effective  dispersion 
coefficient  is  to  inject  a  tracer  into  the  system.  The  tracer  concentration 
is  then  measured  downstream,  and  the  dispersion  coefficient  may  be  found 
from  an  analysis  of  the  concentration  data.  The  first  and  second  moments 
are  usually  needed  to  characterize  a  tracer  distribution  curve.  The  first 
moment  about  the  origin,  which  locates  the  center  of  gravity  of  the  tracer 
curve  with  respect  to  the  origin,  is  usually  called  the  mean  of  a  distribu- 
tion curve  and  the  second  moment  about  this  mean  which  measures  the  spread 
of  the  curve,  is  called  the  variance  of  a  distribution.  One  can  find  the 
functional  relationship  between  the  mean  and  variance  of  the  tracer  distri- 
bution curve  and  the  effective  dispersion  coefficient.  The  formulas  used  to 

2 

evaluate  the  mean,  )^,  and  variance,  (^  ,  of  a  distribution  are 

p   =  f  0  E(6)  d0  (4-113) 


j 


0 


and 
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2=f  (e-^)2  e(8)  de 


or 


=  f  e2s(B)  de-fpe  bCG)  del2  (^-^) 


ex2 

^0 


Using  a  perfect  delta- function  input,  the  tracer  distribution  evaluated  at 
y)  =  i  has  been  found  to  be 


1 


"_  (1-6)  1  (4-10*0 


—■-r-  exp 

V.   7r  t    J  V  T. 


From  this  equation  the  mean  and  variance  are  found  to  be  (22,  23), 

U  =  i  +  2(-2-) 
VXL 


(4-115) 


^2=2(=E_)+8(JL)2  (^-116> 


Solving  for  E/V  L  gives 


(4-117) 


JL  = -L  [  (8*?  +  l)2  -  l] 

If  the  velocity  distribution  of  the  fluid  passing  through  the  system  is 
known,  the  dispersion  coefficient  can  be  obtained  by  injecting  an  impulse 
of  tracer  into  the  system  and  finding  the  variance  from  the  experimental 

data. 

When  S/V^  is  far  smaller  than  1,  Equations  (115)  and  (116)  approximate 

1**1 
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(A. 


=  2(=£-) 


VXL 


which  are  the  sane  as  Equations  (107)  and  (108). 

Consider  a  system  of  steady  state  flow  in  a  straight  pipe  from  the  point 
of  view  of  their  dimensions.  It  has  been  found  (26)  that  the  diraensionless 
group  E/V d  is  a  function  of  the  Reynolds  number,  the  Schmidt  number  and  a 
relative  roughness  number.  The  roughness  factor  has  been  shown  to  be 
important  only  in  turbulent  flow.  The  functional  relationship  for  laminar 
flow  in  a  pipe  may  therefore  be  expressed  as 

dV  P 


V*4 


(  y,      )t(  f>J)) 


(4-118) 


The  generalized  expression  of  the  effective  dispersion  coefficient  is 
E=£-^  (4-119) 


or 


Vxd   *   >*    ^D 

This  relation  is  applicable  only  when  the  system  satisfies  Taylor's  limiting 
conditions,  that  is, 

2 

i-    >>  0.0682  4- 
vm 

for  fluid  flow  in  a  cylindrical  tube  and 
L  %Z 

y"   >  >   0.101  ■%■ 

m 
for  fluid  flow  through  a  slit.  A  family  of  curves  with  E/V  d  as  a  function 
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of  Reynolds  number  and  Schmidt  number  .    8  parameters  is  shown  in  Figure 
(18).  Equation  (119)  can  also  be  writ    as 

From  this  a  family  of  curves  of  slope  1  and  intercept  K/4-  is  resulted  as 
plotted  in  Figures  (19)  through  (21). 

DISPERSION  OF  MATTER  IN  TURBULENT  FLOW  OF  NON- NEWTONIAN  FLUIDS. 

In  a  later  paper  Taylor  (7)  has  shown  that  the  analogous  problem  of 
dispersion  in  turbulent  flow  can  be  solved  in  the  same  way  as  that  used  for 
laminar  flow.  In  the  case  of  Newtonian  flow,  the  effective  dispersion  coeffi- 
cient E  is  found  to  be  10.1  RV*  or  S  =  7.14  RVjJT  (26)  in  which  f  is  the 
Fanning  friction  factor  and  V*  the  friction  velocity.  The  effective  disper- 
sion coefficient  is  almost  entirely  due  to  the  variation  of  the  time-average 
velocity  with  radius  and  the  effect  of  turbulent  axial  diffusion  is  almost 
negligible  in  this  connection.  Taylor's  calculation  assumed  the  validity  of 
Reynolds '  analogy,  and  the  universal  velocity  distribution.  The  turbulent 
flow  region  is  of  the  greatest  practical  importance  in  the  fact  that  most 
fluid  flow  in  industrial  systems  is  turbulent.  It  is  the  purpose  of  this 
section  to  extend  the  analysis  to  non-Newtonian  turbulent  flow  systems.  The 
theoretical  analysis  for  turbulent  flow  of  non-Newtonian  fluids  through  smooth 
round  tubes  was  performed  for  the  first  time  by  Dodge  and  Metzner  (3).  By 
analysis  and  experiments,  they  showed  that  a  purely  viscous  non-Newtonian 
fluid  provides  a  small  amount  of  drag  reduction.  Various  observations  and 
conclusions  regarding  the  different  trends  in  frictional  drag  were  made  by 
many  investigators.  However,  the  turbulent  velocity  profile  suggested  by 
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Fig.     13.     Correlation     for     the     dispersion     of     fluids 
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flowing      in     pipes. 
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(Re)(Sc)=  d  VD 
Correlation     of    dispersion    coefficient    of 
Bingham     plastic     fluid      in     lamina    flow 

in    a    slit. 
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FiqZO.   Correlation    of    dispersion  .coefficiet    of 
Ostwaid-de-    Waele     fluids     in      laminar 

fiov;     in     a     slit. 
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Fig.  21.    Correlation    of    dispersion     coefficient     of 
Bingham    plastic     fluids     In     Icmincr    flow 

in     c     cylndricc'     'cube. 
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Wells  (24)  is  assumed  in  the  present  analysis,  that  is,  the  velocity  V  at 
radius  r  in  a  pipe  of  diameter  2R  is  given  by 


V. 

V*   k 


-£  =  I  In  (l-p  f  (^  (4-122) 


where 

$  =  r/R 

k  =  mixing  length  constant 

V^=  Friction  velocity,  defined  by  (^q/p)2 

Qy=  constant  of  integration 
It  should  be  noted  that  when  k=  0.4,  Equation  (122)  reduces  to  the  expression 
for  Newtonian  fluids  given  by  Prandtl  (25),  which  is  usually  known  as  the 
"universal  velocity  profile".  Reynolds'  analogy  will  also  be  assumed  to  be 
true,  which  ensures  that  the  transfer  of  matter,  heat  and  momentum  by  turbu- 
lence are  exactly  analogous.  Let  e   be  the  coefficient  of  transfer,  Reynolds' 
analogy  in  the  present  case  may  be  expressed  as 

in  which  X   is  the  shear  stress  at  radius  r  and  m  is  the  rate  of  radial 
transfer  of  matter  of  concentration  C.  Two  different  flow  processes,  one  in 
which  a  fluid  flows  through  a  slit  and  the  other  in  which  a  fluid  flows 
through  a  circular  tube  are  treated  separately. 
1.  Flow  through  a  slit 

The  velocity  distribution  in  the  slit  is 

^=|ln  (1-5)+  \  (4-124) 

where 
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4  -  br/i 

and  Reynolds'  analogy  may  bo  expressed  by  the  equation 

From  Equation  (124)  the  mean  velocity  V  may  be  found  as 


x 


\=   V*  J      Vx  d£ 


=  V.  f  (  i  in  (1-^)  t  C,]  d.4 


1. 
0 

=  V*  (Cx   -  1)  (4-126) 

The  viscous  stress  T  at  ^  is  related  to  Tn  by  the  equation 

r  =  i04-1  (4-127) 

Substituting  T  from  Equation  (127)  into  Equation  (125)  yields 

r  4-1 

From  Equations  (124)  and  (128),  one  has 

6  =  kSV^(l  -O  (4-129) 

Based  on  Reynolds '  analogy  the  equation  for  conservation  of  C  is 

±9.  u  -£-2)  =  V  —  +  —  (4-130) 

Introducing  the  dimensionless  variable  ^  =  j  y/& ]  and  substituting  for  £  and 
V  from  Equations  (129)  and  (124),  Equation  (130)  becomes 
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where  x  is  measured  along  the  slit  from  a  fixed  point.  It  is  convenient  in 
the  present  discussion  to  define  concentration  and  velocity  relative  to  the 
axes  -which  move  with  the  mean  flow.  Relative  to  these  axes,  the  velocity  is 

=  y,[ib  (i-U+  1]  (4-132) 

and  Equation  (131)  becomes 

The  small  radial  variation  in  C  can  therefore  be  calculated  from  the  equation 

where  ^C/dx,  is  considered  to  be  independent  of  £.  C  is  then  of  the  form 
C  =  CXi  4-  C;  (^135) 

where  Cx  is  independent  of  £  but  varies  linearly  with  x.,,  and  C»  is 
independent  of  x-,.  Equation  (13^)  can  then  be  solved  as 


a/   1  /"Mlnd-^t  lld^ 


&~  k2  dx- 

k2  °^l  ^o   £ 
k    1  n=l  n 
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The  rate  of  transfer  of  C  across  the  section  at  x,  is 

=  2w4-  Hi  (Si)/1!  J  m  (i-4)+  J]  I   A4n)  d< 

k2    dxi     Jo[k  k  J  n-1    n2 

=    2W*Hf  fe)  f1  ( m  (1-0    f    (-^0+    r     (A4n)]d4 
k3    ^1    ;0    l  n-1    n2  n-1    n2        J 

=    -0.S0S  VJ6—  fe-^)  (4-137) 

k^     ^1 

The  resulting  expression  shows  that  an  effective  coefficient  of  diffusion 
S  would  transfer  matter  across  a  section  at  rate  -2WSE(dC  /cbc),  so  that 
matter  is  transferred  relative  to  a  plane  which  moves  with  its  mean  velocity 
exactly  as  though  it  were  dispersed  by  an  apparent  coefficient  of  diffusion 

S  =  0.404  5Xi  (4-138) 

k3 

In  order  to  estimate  the  effect  of  longitudinal  diffusion  in  turbulent 

diffusion,  it  is  probably  not  unreasonable  to  assume  that  the  coefficient 

of  longitudinal  diffusion  is  equal  to  e.     The  rate  of  transfer  of  matter 

across  a  plane  due  to  longitudinal  diffusion  is 

Q'  -  ZYiS^fi.   f    6  dl  (4-139) 

dx  J0 

Since 

e-kSV^fr-$>  (4-129) 

one  has 

Q»  =  0.333  WS2  kVA(12)  ,  (4-140) 

v   dx 
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which  implies  that  the  mean  coefficient  of  diffusion  due  to  the  longitudinal 
components  of  turbulent  velocity  is 

E1  =  0.16?  k&V^  (4-141) 

Therefore,  the  corrected  value  of  S  which  including  the  longitudinal  turbulent 
diffusion  is  obtained  by  adding  Equations  (138)  and  (141),  that  is, 

E  -  (Mg£  +  0.167  k)  b  V„  (4-142) 

k3 

By  definition 

while  the  Fanning  friction  factor  f  is  defined  by 

r  =  I   f  (O  v2  (4-144) 

0   2   '   x  v     ' 

Thus,  the  friction  velocity  V^  and  the  Fanning  friction  factor  f  can  be 
correlated  by  (26) 


{*.=  /§*  (W*5) 


Therefore,  one  obtains 


E  -  (fiaSgft  -t-  0.167  k)  6  Vv  [if  (4*146) 

^3  x  N  2 

_£_  =  Jl  (Moi  +  0#l6?  k)  fJ-F  (4-147) 

4&VV    a     k3 

■A.  t 

When  k  =  0.4,  that  is,  for  the  special  case  of  the  Newtonian  system  (2) 
the  final  expressions  are 
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E=6.4S&VX  (4-148) 


E—  =1.14  ^T  (4-149) 


Wx 


2.  Flow  through  Cylindrical  Tubes 

The  generalized  velocity  distribution  and  Reynolds '  analogy  have  been 
given  as  (7) 

V    i 

-2E  =  ±  In  (1-|)  4-  C,  (4-122) 

where 
and 

From  Equations  (122)  and  (123)  we  have 

6=kRV^(l-^)  (4-150) 

The  mean  velocity  may  also  be  found  as 

-  V*  (Cx  -  JL)  (4-151) 

By  using  Reynolds '  analogy,  the  equation  for  conservation  of  C  can  be 
written  in  cylindrical  coordinates  as 

-I2  ^r  4s)  =  r<Vx  4s  +  11)  (^152) 

3r     3r      x  9x   et 
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Introducing  dimensionless  variables  and  substituting  for  <£  and  V  from 
Equations  (150)  and  (122),  Equation  (152)  becomes 


3 


jfK<M>  |f  )  =  |^  (l-D  +  <a]K-§|  +  £.§§}  PM53) 


It  is  convenient  to  use  axes  which  move  with  the  mean  speed  of  flow,  which 
is  defined  by 

x1  «■  x  -  Vxt  (4-154) 

Relative  to  these  axes,  the  velocity  distribution  is 

and  the  dispersion  equation  becomes 


The  radial  variation  in  C  can  be  calculated  by  setting  ^  c/ ^  t^  =  0  and 
SC/^x.  being  independent  of  x-,  and  £.  C  is  then  of  the  form 

c  =  0^+  c^  (4-15?) 

where  C  •,  is  independent  of  \   but  varies  linearly  with  x  ,  and  C  c  is 
independent  of  r.  Equations  (156)  is  then  solved  as 


2k2dx. 


V 

^  f  In  (1-P  -  I  In  (1-p  -  1  -  Z    (1  fcn)l     (4-158) 
1  l         I  n=l  n*    J 
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The  rate  of  transfer  of  C  across  the  section  at  x^  is  then 

-i-  Z  c-i,  i))|d| 

n=l  n    J 

3.  1  dcxl 


=  0.323  TjR^  -4* 


clx1! 


(4-159) 


The  effect  of  the  longitudinal  component  of  turbulent  diffusion  can 
be  estimated  by 

•R 

'0 


Q»-  -  2iT(^)  J 


=  -  2*B?V,k  (|j)/  f(l-|)  d| 


=  -  0.16?^  R3kV*  ^S 


dx 


(4-160) 


Q  and  Q'  are  additive.  Comparing  Equations  (159)  and  (160)  with 

Q  =^RTI  -_2^,  the  corrected  value  of  E  allowing  for  longitudinal  diffusion 


is 


ax]_ 


k3 


(4-161) 


Since  the  friction  velocity  V^  and  the  Fanning  friction  factor  f  are 
correlated  by 


V,=  vx 


one  obtains 


V 


(4-145) 


E  „  j£l   (0^28  +  0>l6?  k)  -  j-j 
2    k3  x 


(4-162) 


m 


and 

JL  =  Jfe£I  (9.-328      0#l6?  k)   ;^y  (4-163) 

VXD   *    k3 

For  the  special  case  of  the  Newtonian  flow,  these  expressions  reduce  to 

E  =  5.19  RV#  v  (4-164) 


and 


E 


=  1.84  ,/~f"  (4-165) 


Comparing  Equation  (165)  with  Taylor's  result,  E/V __d  =  3.57 rJ~~ f»  it  is 
seen  that  a  large  difference  exists.  But  if  one  sets  k  =  0.32,  one  can 
obtain  from  Equations  (l6l)  and  (163) 

E  =  10.054  (RV#)  (4-166) 

and 

JL  =  3.56,JT  (4-16?) 

which  are  identical  to  Taylor's  results.  The  velocity  distribution  for 
k  =  0.32,  k  =  0.4  (generally  accepted  for  turbulent  Newtonian  flow)  and 
the  experimentally  measured  velocity  profile  used  by  Taylor  are  plotted  on 
Figure  (22)  for  comparison.  One  conclusion  which  may  be  drawn  is  that  the 
magnitude  of  the  turbulent  dispersion  coefficient  depends  mainly  on  the 
velocity  distribution  near  the  tube  wall.  As  shown  in  Figure  (22)  the 
velocity  profile  used  by  Taylor  and  the  universal  velocity  profile  for 
k  =  0.32  coincide  for  \   in  the  interval  0.9  <$  ^l.  One  can  also  see  from 
the  figure  that  a  large  difference  exists  between  Taylor's  experimentally 
measured  velocity  distribution  and  the  theoretical  result  based  on  Prandtl's 
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mixing  length  theory  (25).  From  the  stand  point  of  fluid  dispersion,  however, 
the  universal  velocity  profile  with  k  =  0.32  is  more  appropriate  in  inter- 
preting the  result  than  the  generally  accepted  value  of  k  =  0.4. 

Levenspiel  (26)  has  mentioned  that  the  molecular  diffusion  in  turbulent 
flow  does  not  contribute  significantly  to  the  transport  of  material.  The 
molecular  diffusion  is  greatly  overshadowed  by  turbulent  eddy  mixing.  Thus, 
a  dimensionless  group  involving  the  longitudinal  dispersion  coefficient  may 
be  expressed  in  terms  of  the  Reynolds'  number  and  the  relative  roughness 
factor,  that  is, 


E 


dVvP 


Vxd"   l  *   '  d 


(4-168) 


In  a  theoretical  analysis  of  mixing  in  turbulent  flow,  E/V  d  has  been  shown 
to  be  a  function  of  the  Fanning  friction  which  is  based  on  the  assumption  that 
Reynolds'  analogy  holds.  In  turn,  however,  f  has  been  experimentally  found  to 
be  a  function  of  the  Reynolds  number  and  the  relative  roughness  parameter. 
Thus,  one  can  see  that  Equation  (166)  is  valid.  A  family  of  curves  which 
correlate  E/vLd,  the  longitudinal  dispersion  number,  and  the  Reynolds  number 
are  shown  in  Figure  (18).  The  curves  for  smooth  tubes  only  and  thus  the 
roughness  factor  is  not  indicated  in  the  figure. 


87 


V.  EXTENSION  OF  ARIS'  WORK  TO  NON-NKWTONIAN  FLUID  FLOW 

If  a  pulse  of  a  tracer  is  injected  into  a  flowing  stream,  this  dis- 
continuity spreads  out  as  it  moves  with  the  fluid  past  a  downstream  measuring 
point.  Levenspiel  and  Smith  (27)  have  shown  that  the  variance  of  a  tracer 
curve  can  conveniently  relate  this  spread  to  the  dispersion  coefficient. 
Taylor's  analysis  has  been  proved  to  be  an  easy  way  to  relate  the  variance  of 
the  tracer  distribution  and  the  effective  dispersion  coefficient.  Since  the 
partial  differential  equation  for  diffusion  cannot  be  solved  in  general,  Aris 
(9)  presented  a  new  approach  to  the  analysis  of  the  problem,  which  calculates 
the  moments  of  the  distribution  in  terms  of  their  respective  parameters,  and 
then  compares  the  moments  with  the  moments  of  the  axial-dispersed  slug  flow 
model  to  find  the  relationship  between  parameters.  His  analysis  removed  the 
restrictions  Taylor  imposed  on  some  of  the  parameters.  The  result  shows  that 
the  rate  of  growth  of  the  variance  is  proportional  to  the  sum  of  the  molecular 
diffusion  coefficient  D  and  the  effective  dispersion  coefficient  suggested  by 
Taylor.  It  is  the  purpose  of  this  chapter  to  extend  Aris '  analysis  to  the 
Ostwald-de  Waele  fluid. 

The  equation  of  continuity  of  the  tracer  in  a  system  with  constant 
density  f>    and  molecular  diffusion  coefficient  D  is  generally  expressed  as 

JL£  =  D  V2C  -  V-  V  C 
5t 

It  is  convenient  to  define  concentration  and  velocity  relative  to  axes  which 
move  with  the  mean  speed  of  flow  and  to  write  the  equation  of  continuity  in 
dimensionless  form.   Consider  a  tube  of  circular  cross  section  and  cylindrical 
coordinates  within  the  tube  (see  Fig.  (2)).  Let 

1  =  (x  -  Vxt)/L,     I   =  r/R,      9  =  tV^L 
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The  governing  diffusion  equation  becomes 

9C=     LP    [1     3n   /fe  ec%   |     1     &2C]  ,      D      32C        ,  2         n-ff  >nflx    a  G 

(5-2) 

which  is  to  be  solved  subject  to  the  initial  and  boundary  conditions 

C(f.  t.  (>.   0)  =  C0(*f,  *.  0),  (5-3) 

±C  „  o  ,       at  5  =  1  (54) 

a5 


Let 


Cp($.  0,  e)  =  J"~Vc(7,  £,  0,  6)  dY  (5-5) 


ana 


mp 


r2TT     rl 
(6)  =  C  =  1   f    d0  /   C(l,   0,  9)  £d$  (5-6) 

p   T[R  70       -'O    P 


be  the  pth  moment  of  the  distribution  of  the  solute  as  a  function  of  },  0 

and  8  and  the  pth  moment  of  the  distribution  of  the  solute  in  the  tube. 

Multiplying  Equation  (2)  by  rp  and  integrating  the  resulting  equation  with 

respect  to  ^  from  -<x»to  «*>,  which  is  performed  by  integration  by  parts, 

yield 

2 


d  8        R2?-    I  S    dI  d  5  S2     d  f  J  LVY  n+1 


'x 


-  ss  rx>  vi      (5-7) 


n+l  *     P 
Equations  (3)  and  (^)  are  also  integrated  to  give 

Cp  (J,  0,  0)=Cp0  (J.  0)  (5-8) 
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-t^-0.  at   J-l  (5-9) 

^5 


The  average  pth  moment  of  the  distribution  can  be  obtained  by  stunning 
Equation  (7)  over  a  cross  section  and  then  dividing  by  the  cross-sectional 
area.  The  use  of  Green's  theorem  (28)  and  the  boundary  condition  imposed 
by  Equation  (9)  ensures  that 

because  Green's  theorem  reduces  the  surface  integral  into  a  line  integral 
along  the  wall  of  the  tube  and  Equation  (9)  provides  ^Cn/«3^  =  0  at  the 
tube  wall.  The  average  pth  moment  of  the  distribution  thus  takes  the  form 

For  p  =  0,  Equation  (10)  gives 

dMQ 

If   =  0  (5-11) 

or  mQ  cs  a  constant,  which  may  be  taken  as  1,  since  the  total  quantity  of  the 
tracer  is  constant.  For  p  =  1,  Equation  (10)  becomes 

dM1   W2TT   r1  o     ^o  n+1 

d0   tt-'q     -/q    n+l   n+1  *     u  3 

For  p  =  2,  Equation  (10)  gives 

A  general  solution  of  Cq  and  C-j_  must  be  obtained  in  order  that  the 
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first  and  second  moments  of  the  tracer  distribution  can  be  evaluated.  Thus 
one  needs  to  solve  Equation  (7)  for  p=  0  and  p  =  1.  For  p  =  0,  Equations  (7) 
and  (3)  reduce  to 

and 

CQ(S.  0.  0)=  C00(J,  0)  (5-15) 

The  solution  of  Equation  (14),  subject  to  the  boundary  condition  ^Cq/^^  =0 
at  \  —  1,  has  been  found  to  be  (9t  20) 

C0(|,  0,  9)=  i.  1+  JL  Z  I  V  cos  m0  +  V  sin  m0] 
K  I    m-0  n=l  u  J 

where  k  =•  LD/R  vx  and  ^^^  is  the  nth  root  of  dJm(^)/d^  =  0.  The  constants, 

A  and  B,  can  be  determined  by  Equation  (15)  and  making  use  of  the  orthogonality 

of  the  functions,  that  is  (28) 

0  2 

/-2TT  yl      Amn    COS   m0    r  -\2 

0  y0      B^  sin  m0  I  J 


r2lj  rl     cos  mj) 

I       dM  ,     Jm^P  cos  <S.0>Sd$  (5-17) 


Substituting  Equation  (16)  into  Equation  (12)  and  integrating  over  the  cross 
section,  one  obtains 

All  the  terms  with  m  =£  0  vanish  by  integration  around  the  circle,  because 
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the  integration  of  a  sine  function  or  a  cosine  function  around  a  circle  is 
zero.  If  n  is  a  positive  integer,  the  integration  can  easily  be  performed. 
Integrating  Equation  (18)  with  respect  to  0  from  0  to  6,  one  obtains 

Kl-  i  Z  4n«-2  f1  -  ■»  (-^2e)]^  2^  -  *2  fen+1)a0<«*|)d|  (5-19) 


As  0—-  0,  }Ln  =  0.  This  means  that  the  center  of  gravity  of  the  tracer 
moves  with  the  mean  speed. of  the  flowing  fluid  if  the  solute  concentration 
is  initially  uniform  over  a  cross-sectional  plane. 
For  p  =•  1,  Equation  (7)  becomes 

3C1_    ED    fl     a    ,t.  *\  .     1    *  °1]  .    .  2        £3^+1:  ,         , 

**~7fJ$  n   il    n]       f    502j+^+l-n+ll       )C0       (5-20) 

in  which  CQ  is  the  same  as  given  in  Equation  (16).  The  complete  solution  of 
C-,  consists  of  two  parts,  one  is  the  complementary  function  which  is  (9»  20) 


e»^>     BO    j»    _  _ 

U  X  £   V  cos  m0  +  B^  sin  m0  j  J^^)  exp  (-ko^B) 
K  m=0  n»l  J 

The  other  is  the  particular  solution  which  may  be  written  as 

oe    cxj 

Substituting  this  equation  into  Equation  (20),   it  is  found  that  *)*„(£)  must 
satisfy  the  equation 


S  H  *£  T  JmWSAn+1       n+1  I 


) 
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l^nen  0  is  large,  the  expression  for  C-,  reduces  to 

where  k  =  LD/R  v  and  K  are  constants.  Substituting  the  value  of  C-j_  into 
Equation  (13)  gives 

dM,   9n    2R27   rl    9    ^^  n+1   2    „   n+3 

_2  =  2P x/    (_2__n+2v   )(k  __L.fc   )*d* 

de   Ly    (n+l)LDJ0    n+1   n+1  3        n+3  5    S 


0  [exp  (-kA)l 


R 


*?, 


=  ^i.  x  __^£  .  , 1 -  (5-22) 

1  10    (n+3)(n+5) 


Lvx 


Thus, 


K2  = 


f  2D_.  +  £  .  1 ]  8  +  0(6)  j-   a  constant        (5-23) 

LVx    LD    (n+2)(n+5)J 


and  the  constant  is  negligible  by  comparison. 

For  the  axial-dispersed  plug- flow  model,  Bischoff  and  Levenspiel 
obtained  the  following  expressions  for  the  moments: 


no- 

1 

*L- 

0 

M2  = 

29 

PL 

(5-24) 

where  l/Pr  is  the  dispersion  coefficient  for  the  axial-dispersed  plug- flow 
model  (17).  Thus  if  V  is  the  variance  of  the  distribution  of  solute  about 
the  moving  origin,  it  is  reasonable  to  define  the  effective  dispersion 
coefficient  as 

S  =  lim  I  %L  (5-25) 

t  —  ~  2  dt 
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It  should  be  noted  that  M  is  the  dimensionless  form  of  V.  By  comparing 
Equation  (23)  with  Equation  (2^),  E  can  be  found  to  be 

R  V 
S=D  + i - — 2S  (5-26) 

2(n+3)(n+5)   D 
which  is  the  sum  of  the  molecular  diffusion  coefficient  D  and  the  apparent 
dispersion  coefficient  suggested  by  Taylor  in  his  first  paper  (6).  But 
there  is  no  restriction  on  the  value  of  their  parameters. 

For  the  case  of  fluid  flow  through  a  slit  with  the  velocity  Vx  everywhere 
in  the  direction  Ox  and  being  a  function  of  y  (see  Fig.  (13) )»  the  equation 
governing  C  is 

The  moving  coordinates  which  move  with  the  mean  speed  of  flow  are  introduced 
in  the  equation  above.  Defining  the  dimensionless  parameters 

*l  =  x/l  ,        i  =  |y/s|.        £  =  s/$,       e  =  tVx/L 

the  governing  diffusion  equation  becomes 


£  C    LD 


58  "7V1  *>2 


SVX 


*>2C  ,  sV  -"»  n+1 


+  ^r    rTr        aM2    W   XlH  *    }  ^         (5"28) 


<H    ar  J  lvx  a>7 


The  constant  of  integration  is  to  be  evaluated  by  using  the  initial  and 

boundary  conditions 


c(»[.  I,  t,   0)-  c0(T.  4.  I)  (5-29) 

"ff  *  °  »  at  {  -  1  (5-30) 

The  pth  moment  of  the  distribution  of  solute  in  the  filament  at  ^,  T  and 


9^ 


at  time  9  and  the  pth  moment  of  the  solute  in  the  tube  are  defined  by 

C  U,  T,  6)  =  /  f?pC(7,  4.  -C>   ©>  d^  (5-3D 


and 


(9)=  C  =  J    dX  /    CL(4.  T,  9)  d£  (5-32) 


24 

Multiplying  Equation  (28)  by  ^p  and  integrating  the  resulting  equation  with 
respect  to  *7  from  -»to<«  yield 

4£  -  "f-  f  A  +  A I  +  ^"1)CP-2  4-  +  P  W  4r  -  SB  ^) 
S0    *%  I  :S| 2   aT 2  -1        ^  LVX    P  X  n+l   n+1  * 

(5-33) 

The  initial  and  boundary  conditions  become 

Gp(4.  T»  o)-  Cp0(5.  X)  (5-3*0 

-jjf  =  0,        at   4  =  1  (5-35) 

The  pth  moment  of  the  solute  in  the  slit  is  obtained  by  averaging  Equation 
(33)  over  the  cross  section.  The  use  of  Green's  theorem  and  Equation  (35) 
insure  that 


0     JQ    c 


1>\        ^2Cn 


n2    ^t2 

Equation  (33)  is  thus  reduced  to 


d^=  0 


de  =  p(p-D  >w  b  +  p j1  dT  r1  (_i_ .  .a  ^  Vi  d^  (5.36) 


'0  J0      n+l      n+l 
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For  p  =  0  Equation  (36)  becomes 


dM 


d0 


°=0 


(5-37) 


or  Mq  =  1  since  the  total  quantity  of  tracer  may  be  taken  as  1.  For  p  =  1 
and  p=  2,  Equation  (36)  gives 


dM. 


d6    -h     ^A    n+l   n+!       ° 


(5-38) 


and 


d8   tv    J0  J0      nvl   n+1       X 


LVV   ^0 


(5-39) 


which  require  general  solutions  of  CQ  and  C-,. 

For  p  =  0  Equations  (33)  and  (34)  reduce  to 


ID 
ft 

X 

^4 

^o  1. 
ar2j 

4 

^2c0 

2  + 

^4 

3T2J 

and 


(5-40) 


c0(4,  *,  o)  =  cQ0U,  r) 


(5-41) 


The  function  satisfying  Equation  (40)  and  given  boundary  conditions  has  been 
given  as  (9»  20) 


c0(4.  t,  e)  =  1  [n-  £  Anvn(4.  t)  exp  (-k\,e)J 


(5-^2) 


where  A  is  a  sequence  of  eigenvalues  and  V"n  is  a  complete  set  of  orthogonal 
eigenfunctions  satisfying  the  equation 
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a  2v„   3  2v. 


St §+\,V-0  (5-43) 


.  2  '     2   n  n 

and  the  boundary  condition 


-£  =  0  at     4-1  (5-^4) 


Inserting  the  value  CQ  into  Equation  (38)  and  integrating  the  resulting 
equation  with  respect  to  0  from  0  =  0  to  0 ,  it  gives 

^=^riVn(i-exp(^na))/o  dr^  (^.nif4  ) 

Vn(4.  T)  d4        (5-45) 
For  p  =  2  Equation  (33)  is 

£C,         T„     r  a2CL         ^2Cn  s  -  0    ntl 


de         S2V    l   2LZ         dT2J-      Vn+1      n+l^        '     0  ^     ' 

If  we  let  LD/$    v x  =  k,  the  complementary  solution  of  Equation  (46)  is 


I  E   A^vn(4,  r)  exP  (-kAne) 


n-1 

in  which  V  is  a  set  of  eigenf unctions  satisfying  Equations  (43)  and  (44), 
The  particular  solution  may  be  written  as 


K    2(ntl)     n+3        n=l 

Substituting  this  equation  into  Equation  (46),  it  is  found  that  tn(^»  "O 
must  satisfy  the  equation 
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.2,,,  .2, 


9t       ***!  n+2>n+l, 


^t-^+^-.jJt.^^Vt.x)  (5-W 


When  0  is  large,  the  expression  for  C,  may  approximately  be  reduced  to 

C,  =   f  f-  -4—r  (1  -  JL.  /n+1)  +  0(exp  (-kA  8)}+  k]  (5-^8) 

1      k  L     2(n+l)  n+3  n     J        J 

where  X  is  a  constant.      Inserting  the  value  of  C-,    into  Equation  (39).   we  have 

dT  =  ^+J0     dTJ0(^"^4        )kl-27nTiy(l-nT3^        } 

+  Ofexp  (-kAn9)jld4 


2- 

-  1D_       *  Vx  4 


LY  ID         3(rn-4)(2n+5) 


Therefore, 

2. 

M    -  f  2D-  +  2  .  & 1 9  +  a  constant 

LLy  LD         3(n+4)(2n+5)J 

Comparing  this  with  Equation  (25),   the  effective  dispersion  coefficient  can 

be  seen  as 

2_ 
i 

2=  D  + 


SVX 


3(n+4)(2n+5)        D 

Aris  (9)  has  shown  that  the  dispersion  model  with  its  smaller  number 
of  parameters  approaches  normality  for  large  time.  Thus  it  is  sufficient  to 
describe  the  solute  distributions  by  considering  the  first  two  moments,  though 
the  approach  to  normality  for  larger  numbers  of  effective  dispersion  coeffi- 
cients is  slow. 
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VI.  OUTLINE  OF  PROPOSED  FUTURE  WORK 

The  treatment  of  solute  dispersion  in  the  previous  chapters  is  confined 
mainly  to  simple  and  generally  useful  tiro- constant  and  three- constant  non- 
Newtonian  models.  Constant  molecular  diffusion  coefficients,  steady  and  iso- 
thermal flow  situations  and  constant  vessel  geometries  are  assumed  in  order 
to  simplify  the  analysis.  On  arriving  at  the  end  of  the  thesis,  the  writer 
is  in  a  position  to  outline  some  of  the  areas  which  are  connected  with  this 
work.  Certainly,  the  methods  used  can  be  extended  to  more  general  problems. 

Mathematical  modeling  is  a  useful  and  time  saving  tool,  but  its  validity 
must  be  verified  by  experiment  since  some  assumptions  have  been  involved  in 
the  derivation.  Radioactive  isotopes,  electrolytes  or  dye  stuffs  may  be  used 
as  tracers.  The  following  factors  should  be  considered  carefully  in  performing 
the  experiment:  roughness  of  pipe,  bends,  elbows  and  valves  in  the  pipe, 
differences  in  physical  properties  between  the  two  mixing  fluids,  and  skewness 
of  the  observed  profile  in  the  analysis  of  experimental  data. 

In  the  derivation  of  the  dispersion  models,  only  the  case  without  solute 
diffusion  across  the  tube  wall  has  been  considered.  Of  course,  this  is  not 
the  only  case  worthy  of  consideration.  Consider  two  phases  flowing  in  co- 
axial annular  regions  with  velocities  and  diffusion  coefficient  varying  only 
with  radial  distance  from  the  common  axis.  The  inner  and  outer  fluids  may  be 
referred  to  as  the  gas  and  liquid  phase  respectively,  for  which  the  resulting 
distribution  of  concentration  can  be  used  to  interpret  certain  gas  adsorption 
and  chromatography  systems. 

Most  investigations  of  fluid  dispersion  have  been  confined  to  simple 
vessel  geometries.  But  many  interesting  problems  arising  in  practical  oper- 
ations are  those  with  non-Newtonian  fluids  in  some  specified  conduits,  such  as 
flow  of  power-law  model  fluids  in  elliptic  tubes  and  fluids  flowing  between 
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cylinders  which  are  not  co-axial. 

It  has  been  pointed  out  (29)  that  Taylor  and  Aris  discussed  only  the 
special  case  of  injecting  the  tracer  into  a  flow  system,  /iris'  more  rigorous 
study  being  restricted  to  a  slug  of  material  with  initial  finite  moments  with 
respect  to  the  axial  coordinate.   It  is  anticipated  that  simpler  and  more 
convenient  formulations  can  be  developed  from  the  properties  of  the  exact 
solutions  of  the  diffusion  equation.  It  is  obvious  that  a  full  solution  of 
the  diffusion  equation  would  be  very  difficult.  'What  we  desire  is  an  under- 
standing of  the  gross  features  of  dispersion.  The  variation  of  C(x,  t) 
frequently  provides  the  information  desired.  Thus,  the  formulation  in  terms 
of  apparent  Fick's  diffusion  down  a  gradient  of  C  relative  to  the  main  motion 
is  quite  convenient  for  interpreting  systems. 

Much  interest  has  been  given  to  the  flow  of  blood  in  the  human  cardio- 
vascular system  (30,  31.  32).  Blood  is  apparently  non-Newtonian  even  though 
it  has  been  considered  to  be  Newtonian  in  most  investigations.  A  study  of 
the  blood  flow  system  may  be  performed  by  considering  the  dispersion  of  non- 
Newtonian  fluids  associated  with  pulsating  flow  in  tubes  in  which  the  walls 
have  appreciable  flexibility.  Kapner  (33)  bas  studied  such  a  system  by  means 
of  dimensional  analysis.   The  change  of  the  pressure  pulse  resulting  from 
changes  in  the  dimensionless  parameters  has  been  studied. 

A  chemical  engineer  specifying  mixing  equipment  needs  to  predict  the 
flow  pattern,  average  residence-time  and  residence-time  distribution  in  case 
of  continuous  operations.  A  discussion  of  the  literature  prior  to  1959  is 
presented  by  Bernhardt  (3*0  and  I-!cKelvey(35)«  Significant  progress  has  been 
made  in  the  ana-lysis  of  laminar  mixing  processes  by  the  concept  of  shear 
deformation.  Such  theoretical  analyses  are  described  by  Spencer  and  V/iley 
(36),  Kohr,  et  al.(37)  and  Schrenk,  et  al.(3S).  An  annular  channel(38)  and 
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a  tube  (39)  "which  continuously  mix  two  polymer  streams  by  rotation  of  the 
boundaries  to  produce  the  shear  deformation  are  considered.  Good  agreement 
was  reported  between  theory  and  experiment  for  the  mixing  pattern.  For  this 
type  of  mixing,  performance  data  are  needed  for  the  various  types  of  mixers 
available  for  mixing  or  blending  Newtonian  and  non-Newtonian  fluids  of  various 
viscosity  ratios  under  laminar  and  turbulent  flow  regimes  for  various  flow 
ratios  of  components. 
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NOMENCLATURE 

C  concentration 

CQ  concentration  of  tracer  in  inlet  stream 

C  area  mean  concentration 
ex 

C  pth  moment  of  the  distribution  of  solute  in  the  filament  through 

\,   0  at  time  Q 

C  pth  moment  of  the  distribution  of  solute  in  the  tube 

D  molecular  diffusivity,  1  /t 

d  inside  diameter  of  round  tube,  1 

E  effective  or  apparent  dispersion  coefficient,  1  /t 

E(8)  exit  age  distribution  function 

F(B)  response  curve  to  a  step  function  input  of  tracer 

f  Fanning  friction  factor 

K  effective  dispersion  coefficient  factor 

k  mixing  length  constant 

k  dimensionless  parameter  LD/R'TX 

L  characteristic  length  of  system,  1 

Mp  pth  moment  of  the  distribution  of  solute  in  the  tube 

m  parameter  in  Ostwald-de  Waele  model 

n  flow-behavior  index  in  Ostwald-de  V,raele  model 

?  P-pgz 

P,  dispersion  coefficient  for  the  axial-dispersed  plug- flow  model 

p  pressure,  m/lt 

Q  volumetric  flow  rate,  1^/t 

R  tube  radius,  1 

r  radial  distance  from  center  of  tube,  1 
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ro  V(po  -  pl>/2l 

t  time,  t 

t  mean  residence-time,  t 

U  a  unit  step  function 

V  center  line  velocity,  l/t 


m 


V^  slip  velocity  at  tube  wall,  l/t 

V  axial  velocity,  l/t 

V  mean  velocity  in  axial  direction,  l/t 
V#  friction  velocity  defined  by  (Tq/^)2 
x  axial  distance,  1 

x..  axial  position  relative  to  a  coordinate  system  moving  with  the 

mean  speed  of  flow,  1 

y  cross-sectional  distance  in  rectangular  coordinates,  1 

^0  V(VPL>L 

Greek  Letters 

d  n-th  eig en- value 

n  to 

c<  diraensionless  constant 

P  dimensionless  constant 

f  wall  slip  velocity  factor 

A  rate  of  deformation  tensor 

5  characteristic  length  in  the  cross-sectional  direction  in 
rectangular  coordinates,  1 

5(8)  Dirac  delta  function  input 

6  coefficient  of  transfer 

£>  dimensionless  position  variable  defined  by  y/$ 

^  dimensionless  flow  behavior  index  in  Bingham  plastic  model 
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Y]  dimensionless  axial  distance,  x/L 

^,  dimensionless  axial  distance  relative  to  a  coordinate  system 

moving  with  the  mean  speed  of  flow 

6  dimensionless  time,  t/t  or  tV_./L 

A  dimensionless  constant 

^  mean  of  the  r.  t.  d.  f.  or  the  first  moment  about  the  origin 

u.  parameter  in  the  Bingham  model,  m/lt 

V  parameter  in  the  Ostwald-de  Waele  model 

\  dimensionless  radial  distance,  r/R 

\~  dimensionless  flow  behavior  index  in  Bingham  plastic  model 

P  fluid  density,  m/l^ 

2 

oc  variance  of  the  r.  t.  d.  f. 

~C  dimensionless  distance,  z/fc 

X  shear  stress  tensor,  m/t  1 

X"0  parameter  in  Bingham  plastic  model,  m/tl 

*$>  angle  between  the  vertical  direction  and  the  x-coordinate 

yo  parameter  in  Ellis  model 

i^-,  parameter  in  Ellis  model 
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A  study  has  been  made  of  the  dispersion  of  non-Newtonian  fluids  in 
continuous  flow  systems.  If  a  pulse  of  tracer  is  injected  into  a  flowing 
stream,  this  discontinuity  spreads  out  as  it  moves  with  the  fluid  past  a 
downstream  measurement  point.  The  amount  of  spreading  depends  on  the 
intensity  of  dispersion  in  the  system.  It  has  been  proved  that  the  variance 
of  the  tracer  distribution  and  the  effective  dispersion  coefficient  can  be 
easily  related  by  means  of  Taylor's  analysis. 

Dispersion  by  convection  alone,  which  is  usually  called  the  convective 
model,  is  considered  first.  Discontinuous  points  do  occur  in  the  resulting 
distribution  curves.  The  deformation  in  the  shape  of  the  actual  residence- 
time  distribution  curve  is  caused  by  the  dispersion  effect. 

The  dispersion  model  assumes  that  the  combined  effects  of  the  variation 
of  axial  velocity  with  cross-sectional  position  and  the  cross-sectional 
material  transport  by  molecular  diffusion  are  the  dominant  factors  causing 
the  dispersion.  Conditions  under  which  such  assumptions  are  valid  are 
investigated.  Expressions  for  the  residence-time  distributions  as  a  cor- 
relation of  the  dispersion  coefficients  are  obtained. 

An  alternate  approach  is  that  which  characterizes  the  dispersion  in 
terms  of  its  moments.  These  moments  are  related  simply  to  measurable 
quantities.  This  approach  is  used  to  interpret  the  dispersion  of  the 
Ostwald-de  Waele  fluid. 


